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Abstract. We consider solutions u to the 3d nonlinear Schrodinger equation idtu+ 
Au + jitpM = 0. In particular, we are interested in finding criteria on the initial 
data uo that predict the asymptotic behavior of u{t), e.g., whether u(t) blows- up 
in finite time, exists globally in time but behaves like a linear solution for large 
times (scatters), or exists globally in time but does not scatter. This question has 
been resolved (at least for data) in [151 [i H [12] if M[u]E[u] < M[Q]E[Q], 
where M[u\ and E[u\ denote the mass and energy of u, and Q denotes the ground 
state solution to —Q + AQ + |QpQ = 0. Here we consider the complementary 
case M[w]£'[m] > A/[(5]i?[(5]. In the first (analytical) part of the paper, we present 
a result due to Lushnikov [5^, based on the virial identity and the uncertainty 
principle, giving a sufficient condition for blow-up. By replacing the uncertainty 
principle in his argument with an interpolation-type inequality, we obtain a new 
blow-up condition that in some cases improves upon Lushnikov's condition. Our 
approach also allows for an adaptation to radial infinite-variance initial data that 
has a conceptual interpretation: for real-valued initial data, if a certain fraction 
of the mass is contained within the ball of radius ^/[u], then blow-up occurs. We 
also show analytically (if one takes the numerically computed value of ||Q||^i/2) 
that there exist Gaussian initial data uq with negative quadratic phase such that 
||wo|Ihi/2 < ||Q|Ihi/2 but the solution u(t) blows-up. In the second (numerical) 
part of the paper, we examine several different classes of initial data - Gaussian, 
super-Gaussian, off-centered Gaussian, and oscillatory Gaussian - and for each class 
give the theoretical predictions for scattering or blow-up provided by the above 
theorems as well as the results of numerical simulation. On the basis of the numerical 
simulations, we formulate several conjectures, among them that for real initial data, 
the quantity ||Q||^i/2 provides the threshold for scattering. 



The nonlinear Schrodinger equation (NLS) or Gross-Pitaevskii equation is 



with wave function u = u{x, t) G C. We consider x G M" in dimensions n = 1,2, or 3. 
The initial-value problem is locally well-posed in (see Cazenave [5] for exposition 
and references therein). In this now standard theory obtained from the Strichartz 
estimates, initial data uq G give rise to a unique solution u{t) G C([0, T]; H^) with 
the time interval [0,T] of existence specified in terms of HMolli/i- In some situations, 
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an a priori bound on can be deduced from conservation laws which imphes 

the solution u{t) exists globally in time. On the other hand, we say that a solution 
u{t) to NLS blows-up in finite time T* provided 

(1.2) hm ||VM(t)|U2 = +oo. 

For = 1, all initial data yield global solutions, but large classes of initial data 
leading to solutions blowing-up in finite time are known for n = 2 and n = 3. NLS 
arises as a model of several physical phenomena. We outline three important examples 
in a supplement to this introduction ( ^1.1| below), and emphasize that in each case 
the mathematical property of blow-up in finite time is realistic and relevant. It is 
therefore of interest to determine mathematical conditions on the initial data that 
guarantee the corresponding solution will blow-up in finite-time and conditions that 
guarantee it will exist globally in time. Moreover, if we know the solution is global, 
it is natural to ask whether we can predict the asymptotic {t +oo) behavior of the 
solution. If the solution asymptotically approaches a solution of the linear equation, 
we say it scatters. Nonlinear effects can persist indefinitely, however; for example, 
leading to formation of solitons or long-range modulation of linear solutions. 

Partial answers to the above mathematical problem are known, and we will discuss 
separately the existing literature in the case of dimensions n = 1, 2, and 3. Afterward, 
we will state our new findings in the n = 3 case. 

Before proceeding, we note that NLS satisfies conservation of mass M[m], momen- 
tum and energy E[v\^ where 

|2 



M[u] = \\u\\l2, P[u] = Im / uVu 



E[u] = ^llVull^^ - ^11^11^4. 



Also, NLS satisfies the scaling symmetry 

(1.3) u{x,t) solves NLS =^ \u{\x,\'^t) solves NLS. 

Consequently, the critical (scale-invariant) Sobolev space if'*(]R"') is s = The 
NLS equation also satisfies the Galilean invariance: For any v G M", 

u{x, t) solves NLS =^ e'''-''e-''\''\\{x - 2vt, t) solves NLS , 

and thus, any solution can be transformed to one for which P[u\ = 0. Let 

V[u]it) = \\xuit)\\l, 

denote the variance. Assuming 1^[m](0) < oo, then the virial identities (Vlasov- 
Petrishchev-Talanov |2H], Zakharov [32], Glassey [IB] ) 

(1.4) dtV[u]=Umjx-Vuudx, d^V[u\ = 8nE[u] + {8 - 4:n)\\Vu\\l2 
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hold. Let Q = Q{x) denote the real- valued, smooth, exponentially decaying ground 
state solution to 

(1.5) -Q + AQ + Q^ = 0. 

Then u{x,t) = e**Q(x) solves NLS, and is called the ground state soliton. The Po- 
hozhaev identities are 

n 4 

(1-6) IIVQili. = ^\Q\\h , WQWl = ^WQWh ■ 



Weinstein [30] proved that the Gagliardo-Nirenberg inequality 

iri|v0||^ 



(1-7) \m% < CGNii0iiir IIV0II 



is saturated hj cf) = Q, i.e.. 



CGN - „^,|4-n| 



L4 



IIQIIi^iivQii^ 

is the sharp constant. 

2d case. Much of the mathematically rigorous literature has been devoted to the 



2d case, of particular relevance to the optics model (item 1 in [:1.1), and has the 



special mathematical property of being L^-critical. The energy E[u] conservation 



combined with the Weinstein inequality (1.7) implies that if ||mo||l2 < ||Q||l2, an 
solution is global. This result is in fact sharp, in the following sense. The scale- 
invariance of the 2d equation allows for an additional symmetry, the pseudo-conformal 
transformation 

(1.8) m(x, t) solves 2d-NLS =^ m(x, t) = ^ eT^M(^|, ^) solves 2d-NLS. 
This gives rise to an explicit family of blow-up solutions 



X 



(T-t) ^ \T-t 

obtained by the pseudoconformal transformation, time translation, and scaling. They 
blow-up at the origin at time T > (and T can be taken arbitrarily small), but 
\\ut\\l^ = ||Q||l2. Note that they have initial data {ut)o{x) = e*l^l'/^g(a;/T)/T, 
indicating that the inclusion of a quadratic phase prefactor can create finite-time 
blow-up. Moreover, it was observed by Vlasov-Petrishchev-Talanov [28], Zakharov 
[33] and Glassey [13] that if the initial data has finite variance HxmoHlz < oo and 



E[u] < (which implies by (1.7) that ||mo||l2 > ||<5||l2), then the solution u{t) blows- 
up in finite time. Blow-up solutions with E[u] > exist and global solutions with 
E[u] > exist. If E[u] > 0, then a sufficient condition for blow-up can be deduced 



4 JUSTIN HOLMER, RODRIGO PLATTE, AND SVETLANA ROUDENKO 

from the virial identity (see [22], [55] 

(1.9) Vt{0) < -^/l6EV{0). 



Id case. The Id case is subcritical; energy conservation and (1.7) prove that 
solutions never blow-up in finite time. One can still ask if there is a quantitative 
threshold for the formation of solitons. Such a threshold must be expressed in terms 
of a scale-invariant quantity, and the norm is a natural candidate. We note that 
soliton solutions 

u{x,t) = e'^Q{x) , Q{x) = v^sechx 

have 1 1 ^(t) 1 1 = IIQIIli = "s/S^r (as do rescalings and Galilean shifts of this solution). 
As the equation is completely integrable (see Zakharov-Shabat [32]), one has available 
the tools of inverse scattering theory (1ST). 1ST has been applied by Klaus-Shaw [21] 
to show that if HmoIIli < then no solitons form. See Holmer-Marzuola- 

Zworski [16], Apx. B for a calculation showing that this is sharp - for initial data 
uo{x) = aQ{x) with a > ^, a soliton emerges in the t +oo asymptotic resolution. 
We remark that although no solitons appear if ||mo||li < such solutions do 

not scatter, i.e., they do not approach a solution to the linear equation as t — > +oo 
- see Barab [Ij. In fact, there are long-range effects and one conjectures modified 
scattering - see Hayashi-Naumkin [TH] for some results in this direction for small intial 
dataj^ The Hayashi-Naumkin paper, in fact, treats a more general equation and does 
not rely on 1ST; presumably 1ST could be applied to prove modified scattering for 
II'^oIIli < ^IIQIIli for generic Schwartz uq, although we are not aware of a reference. 

3d case. We have previously studied the 3d case of NLS, which is supercritical, 
in Holmer-Roudenko [IZl UHl HH] , Duyckaerts-Holmer-Roudenko [H], and Duyckaerts- 
Roudenko [^. Scattering and blow-up criteria are most naturally expressed in terms 
of scale invariant quantities, and natural candidates are the norm and the H^^^ 
norm. We argue below that the norm is completely inadequate, and while the 
^1/2 norm is a more reasonable choice, it too appears deficient. In [HI [8], we work 
instead with two scale- invariant quantities: M and 

def ||M(i(;) 11^2 II VM(t) 11^2 



(1.10) r]{t) 



\\Q\\l4'^Q\\l^ 



By the Weinstein inequality (1.7) and the Pohozhaev identities (1.6) we have 



2 ^ M[u]E[u] 



^Blow-up solutions are also possible when E — provided Vt(0) < 0, for a general review refer to 



2 



[T5] does not cover the full range ||wo||li < fllQIUi, and, in fact, the smallness condition is in 
terms of a stronger norm. 
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This results in two "forbidden regions" in the M[M]ii^[M]/M[Q]ii^[Q] versus yf phase- 
plane - see the depiction in Figure [Llj Note that since M[u\ and E[u\ are conserved, 
all time evolution in Figure occurs along horizontal lines. In what follows for 
brevity and simplicity we assume P[u\ =0 which can be obtained via Galilean trans- 
form. The most general case would follow as it is explained in Appendix B of |19j . 



Theorem 1.1 (Duyckaerts-Holmer-Roudenko [8], Holmer-Roudenko [HI [19]). Sup- 
pose that uo e and M[u]E[u] < M[Q]E[Q]. 

(1) If ri{0) < 1, then u{t) is globally well-posed and, in fact, scatters in both time 
directions. 

(2) lfr]{0) > 1 and either uq has finite variance oruo is radial, then u{t) blows-up 
in finite positive time and finite negative time. 

(3) // ri{0) > 1, then either u{t) blows-up in finite forward time or there exists 
a sequence tn +cxd such that || VM(t„,)||L2 = oo. A similar statement for 
negative time holds. 

It is a straightforward consequence of the linear decay estimate that scattering 
solutions satisfy 

lim ||u(t)||LP = 0, 2 < p < 6. 

t/'+oo 



It follows by using the p = 4 case and the Pohozhaev identities (1.6) that 

lim .it? = ^M^M 
iiTc«^^^ 3M[Q]E[Q] • 

That is, in Figure |l.l| a scattering solution has ri{t? asymptotically approaching 



boundary line ABC. On the other hand, since blow-up solutions satisfy (1.2), such 



solutions go off to right (along a horizontal line) in Figure |1.1[ We note that Merle- 



Raphael [21] strengthened (1.2): they proved that if u{t) blows-up in finite forward 
time T* > 0, then 

lim \\u(t) \\r3 = -\-oo . 

t/T* 

The blow up for finite variance as in Theorem |1.1[ part (2), has previously been 
obtained by Kuznetsov et al. in |23j . 

The results of Duyckaerts-Roudenko [H] are contained in the next two theorems. 
First, they establish the existence of special solutions (besides e^^Q) at the critical 
mass-energy threshold. 

Theorem 1.2 (Duyckaerts-Roudenko |9j). There exist two radial solutions and 
Q~ of NLS with initial conditions such that Qq E risyoH'^CR.^) and 

(1) M[g+] = M[Q-] = M[Q], E[Q+] = E[Q-] = E[Q], [0, +oo) %s m the (time) 
domain of definition of and there exists Cq > such that 

Vt>0, \\Q^{t)-e''Q\\^,<Ce-'^\ 
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Figure 1.1. A plot of M[u\E[u\/M[Q]E[Q] versus ^f, where r] is de- 



fined by (1.10). The area to the left of line ABC and inside region ADF 



are excluded by ( |1.11 ). The region inside ABD corresponds to case (1) 



of Theorem |1.1| (solutions scatter). The region EDF corresponds to 
case (2) of Theorem |1.1| (solutions blow-up in finite time). Behavior 
of solutions on the dotted line (mass-energy threshold line) is given by 
Theorem 11.31 

(2) llVQo II2 < IIVQII2, Q~ is globally defined and scatters for negative time, 

(3) ||V(5([||2 > llVQIb; CLnd the negative time of existence of is finite. 

Next, they characterize all solutions at the critical mass-energy level as follows: 

Theorem 1.3 (Duyckaerts-Roudenko [9]). Let u he a solution of NLS satisfying 
M[u]E[u] = M[Q]E[Q]. 

(1) If TjiQ) < 1, then either u scatters or u = Q~ up to the symmetries. 

(2) If rjiOi) = 1, then u = e^^Q up to the symmetries. 

(3) 7/77(0) > 1, and uq is radial or of finite variance, then either the interval of 
existence of u is of finite length or u = up to the symmetries. 
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A recent result of Beceanu [2| on ( 1.1 ) states that near the ground state sohton (and 
its Gahlean, scahng and phase transformations) there exists a real analytic (center- 
stable) manifold in H^/'^ such that any initial data taken from it will produce a global 
in time solution decoupling into a moving soliton and a dispersive term scattering in 

In part I of this paper, we provide some alternate criteria for blow-up in the spirit 
of Lushnikov [20]. First, we state his result, adapted to our notation. Due to the 
complexity of the formulas, we will write M = M[u], E = E[u], etc. For simplicity we 
restrict to the case E > 0, since E < is comparately well understood from Theorem 



Theorem 1.4 (adapted from Lushnikov [20]). Suppose that uq G and HxmoIIlz < 
oo. The following is a sufficient condition for blow-up in finite time: 



(1.12) 

where 



M 



3M2 J 



(1.13) 




+ 00 



ifO <LJ <1 



^ +LJ -3 ifuj>l, 



which is graphed in Figure 1.2 



For an explicit formulation of the condition (1.12) refer to ^4.3, in particular, when 



condition rewrites as in (4.6)-(4.7). 



initial datum is real- valued (1.12) becomes (4.4) and when it is complex- valued the 



Theorem 1.4 is based upon use of the uncertainty principle. 



:i.i4) 



9 



Im I ix ■ Vu)u dx 



< 



4 
9' 



xu\\\2 1 



L2 5 



the virial identity, and a "mechanical analysis" of the resulting second-order ODE in 



V{t). By replacing (1.14) with 



:i.i5) 



22.75 



35-52 



3 



4 

I 



we can obtain a different condition which in some cases improves upon Theorem 



1.4 The inequality (1.15) can be thought of as a variant of the Holder interpolation 
inequality \\u\\l2 < \\u\f/J^r^\\u\\t^J , since ||M||i6/5 and ||xM||i2 scale the same way. Both 



inequalities (1.14) and (1.15) are stated here with sharp constants and are proved in 

m 
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Figure 1.2. A plot of g{uj) versus a;, where g is defined in ( 1.13 ). This 



function appears in the blow-up conditions in Theorems 1.4 and 1.5 



Theorem 1.5. Suppose that uq G and ||a;Mo||L2 < oo. The following is a sufficient 
condition for blow-up in finite time: 



(1.16) 



Vt{0) 2V2{ME)h f^CfEl^.. 
< ^ 9 I — 77^^(0) 



M 



Ml 



C 



2^ ■ 7^ ■ TT- 

35-52 



where g is defined in (1.13) and graphed in Figure 1.2 



For an explicit reformulation of (1.16) refer to ^4.3, in particular, for the real- valued 



initial datum it becomes (4.5) and for the complex-valued datum it is equivalent to 



(4.12) - (4.13). 



Note that (1.16) can be put into the form 
'SEV{fd)\ 



M 



< 



2^3, 



3M2 



g{uj)=^g{fi uj) 



V2{ME)^ 



which offers a comparison between Theorem 1.4 and 1.5 These conditions should be 



compared to the sufficient condition for E[u] > in the 2d case, namely (1.9). 



The approach via the interpolation inequality (1.15) also allows us to prove a 



radial, infinite-variance version of Theorem [1.5 based upon a local virial identity, 
Strauss' radial Gagliardo-Nirenberg inequahty [25], and a bootstrap argument. Select 
a smooth radial (nonstrictly) increasing function iIj{x) such that ip{x) = for 
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< < 1 and = 2 for |x| > 2. Define the localized variance 



;i.i7) 



V, 



def 



R 



R^i){x/R)\u{x)\^ dx. 



Note that by the dominated convergence theorem, for any uq & L"^, we have 

VRiO) 



hm 



0. 



Thus, there always exists R such that (1.18) below holds. 

Theorem 1.6. Suppose ME > 1. Fix 6 <^ 1 (the smallness depends only on iIj{x) in 
[1.17]). Given Uq G radial, take any R such that 

Vr{0) 



(1.18) 



M 



R'> 



6 



(the implicit constant in the second inequality again depends only on ijj{x) in {1.17)). 
Then the following is a sufficient condition for blow-up in finite time: 



+ 5)5(1 -5)5(ME)6 {8 + 6)HCoo)^ eI 



(1.19) 



M 



(a 



9 



(1-5) 



Mi 



where g is defined in (1.13) and graphed in Figure 1.2 



One way to generate examples of uq satisfying the hypotheses of Theorem 1.6 but 



not Theorem 1.5 is to take any of the examples detailed below for which Theorem 1.5 



applies, but tack on a slowly decaying tail of infinite variance but at very large radii. 
For example, redefine mq at very large radii to be uo{x) = \x\~'^. However, the main 
merit of Theorem 1.6 is the availability, in the case of real initial data, of a conceptual 
interpretation in terms of the way in which mass is initially distributed. 

Corollary 1.7. There zs < 5 ^ 1 such that the following holds. Suppose that 
ME > 1, Mo G is radial and real, 

1 



(1.20) 



M 



|x|>5i/2M(M£;)-i/3 



\uo\'dx < 6\ME) 



-2/3 



Then blow-up occurs in finite time. 



Note that the quantity on the left-side of (1.20) is the fraction of initial mass 
occurring outside the ball of radius 5M(ME)~^^^. Thus, (1.20) states that most 
mass is inside the ball of radius SM{ME)~^^^, and we intuitively expect an initially 
highly concentrated real solution to blow-up. By the scaling (1.3), it is natural that 
the radius scales linearly with M. 
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Theorems |1.4[ |1.5[ and |1.6[ and Corollary |1.7 are proved in 

is in 



and the reformula- 



tion of the conditions of Theorems |1.4| and |1.5 
In the remainder of the paper 



B]-[9} we examine several specific radial initial data 
given as profiles with several parameters]^ In each case, we report the predictions 
given by Theorems |1.4 and L5, and also the results of numerical simulations. 
In general, Theorems |1.4| [l75 may not give better results than Theorem (we have 
one example where it does not, see Figure 5.3), however, we show that they give new 



information in many other cases, in particular, when initial data has a negative value 
of \4(0). 

In ^ we consider initial data 

where Q is the ground state solution to ( 1.5 ). The case 7 = is completely understood 
by Theorem In fact, \^^'^Q{Xr) corresponds to the parabolic-like boundary curve 
in Figure [lT] which lies below the M[M]ii^[M] = M[(5]-E[Q] line for all A 7^ 1, and we 
have blow-up for A > 1 and scattering for A < 1. The case 7 7^ is more interesting; 



in particular, for negative phase, 7 < 0, Theorems 1.4 and 1.5 give new range on 
parameters (A, 7) for which there will be a blow up, see Figure 5.2 Although for 



positive phase, 7 > 0, there is no new information on blow up other than provided 



by Theorem lA, we show the numerical results for the blow up threshold in Figure 
5.3[ in particular, there is no blow up for A < 1 is expected. 
In ^ we consider Gaussian initial data 



Uo[X) 



pe 



By scaling, it suffices to consider 7 = 0, In the real case 7 = 0, the behavior will 
be a function of p/ 
appears in Figure 



6.1 



6.2 



a, and the results are depicted in Figure 
and the case 1 = —\ appears in Figure |6.3 



The case 7 = | 
For this initial 



data Theorem 1.4 gives the best range for blow up, although Theorem 1.5 gives an 



improvement over the Theorem 1.1 
In ~ 



we consider "super-Gaussian" initial data 

Ui^{x) = pe 



7.1 



By scaling, again it suffices to consider 7 = 0, ±2. In the real case 7 
will be a function of p/ a^^^, and is depicted in Figure 
in Figure 7.3, and the case 7 = — ^ is presented in 
Theorem L5 gives the best theoretical range for blow up 
In ® we consider "off-centered Gaussian" initial data 



r igure 



0, the behavior 
The case 7 = | is presented 
For this initial data 



7.2 



pr'^e 



9 ■ 9 



■^Since we work exclusively with radial data, we write our functions as functions of r €E (0, 
but keep in mind that we are studying the 3d NLS equation. 



-00 
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By scaling, again it suffices to consider 7 = 0, In the real case 7 = 
will be a function of "pjo?!^^ and the results are presented in Figure 

and the case 7 = • 



7 



+ 2 is given in Figure 



8.2 



8.3 



— 2 is given in Figure 
initial data Theorem 1.5| gives as well the best theoretical range for blow up. 
In ^ we consider "oscillatory Gaussian" initial data 



0, the behavior 
The case 
For this 



Uo{x) = p cos(/5r) 



9 ■ 9 



9.3 



respec- 



We restrict our attention to 7 = 0, presented in Figures 9.1, 9.2| and 
tively. For the oscillatory Gaussian the best theoretical range on blow up threshold 
is provided by a combination of Theorems 1.4, 1.5 for small oscillations, (3^1, The- 
1.4 is stronger, and for fast oscillations, j3 > 1, Theorem 1.5 provides a better 



orem 

range (for exact values see the above Figures). 

The numerics described in ^|5]-[9] provide evidence to support the following conjec- 
tures. 

Conjecture 1. For each e > 0, there exists radial Schwartz initial data uq for which 
M > M[(5]-E'[Q], II Mo ~ Q II //I < u(t) scatters as t ^ — 00, and u(t) blows-up 
in finite forward time. 



That is, there exist initial data arbitrarily close to D in Figure IJ^ with the 
property that the backward time evolution results in scattering but the forward 
time evolution results in finite time blow-up. The numerical evidence of the ex- 
istence of such solutions is a consequence of the study of initial data of the form 

Take A < 1 but close to 1. Then we find that there exists a 
curve 7o(A) such that limA/i 70(A) = with the following property: If I7I > 70(A), 
then M[M]ii^['u] > M[(5]-E'[<5], and uq evolves to a solution u(t) blowing up in finite 
positive time if 7 < —70(A) but uq evolves to a scattering solution in positive time if 
7 > 7o('^) (see Figures 5.2 and 5.3). By the time reversal property 

u{t) solves NLS =^ solves NLS 

we conclude that if 7 < —70(A), u{t) scatters backward in time. We can take A as 
close to 1 and 7 as close to as we please (while maintaining 7 < —70(A), establishing 
the claimed conjecture (as a result of observed numerical behavior). 

Before proceeding, let us remark on some consequences assuming this conjecture is 

that the smallest admissible value of ?](0)^ that could 



1.1 



valid. We see from Figure 
lead to a finite-time blow-up solution is |-|- (Corner B). 

1st Corollary of Conjecture 1. For each e > 0, there exist initial data Uq with 
M > M[Q]E[Q] and ■r]'^{0) < ^ + e for which the evolution u{t) blows-up in 
finite time. For each N ^ 1, there exists initial data uq with M[-u]i?[M] > M[(5]-E[(5] 
and ?7^(0) > leading to a scattering solution u(t). 



More loosely stated, there exist initial data as close to point B in Figure 1A_ leading 
to finite-time blow-up solutions, and there exist initial data as far to the right in 
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the direction of E in Figure 1.1 leading to scattering solutions. This establishes 
the irrelevance of the size of 77(0) in predicting blow-up or scattering in the case 
M[u]E[u] > M[Q]E[Q]. 

This follows from Conjecture 1 as follows. Taking u{t) to be a solution of the 
type described in Conjecture 1, note that \imt\_oori'^it) = |. Hence, for some large 
negative time — T, we have | < rfi^—T) < | + e. Resetting, by time translation, to 
make time — T into time gives the first type of solution described here. On the other 
hand, if T* denotes the blow-up time of u{t), then \im.t/T* v'^{t) = +00. Therefore, 
there exists a time T' < T* but close to T* such that 1]'^{T') > N. Applying the time 
reversal symmetry and time translation gives the second type of solution described 
here. 

2nd Corollary of Conjecture 1. For each e > 0, there exists initial data Uq with 
II'^oIIl^ < e for which u{t) blows-up in finite time. For each ^ 1, there exists initial 
data Uq for which HmqIIls ^ ^ and for which u{t) scatters. 

Stated more loosely, the quantity ||mo||l3 is irrelevant to predicting blow-up or scat- 
tering. Thus, the critical Lebesgue norm gives no prediction of dynamical behavior 
in the 3d case; note the contrast with the Id case discussed above, where the critical 
Lebesgue norm determines the threshold for soliton formation. 

This follows from Conjecture 1 by the same type of reasoning used to justify the 
first corollary, since if u{t) scatters in negative time we have limf\_oo ||'u(i(:) ||/^3 = 
and if u{t) blows-up in finite forward time T*, we have limtyT* \\u{t)\\L'J' = +00. This 
latter fact was proved by Merle- Raphael \24\ . 

What about the H^^'^ norm? The "small data scattering theory" (essentially a 
consequence of the Strichartz estimates - see for exposition) states that there 
exists S > such that if ||uo||^i/2 < 6, then u{t) scatters in both time directions. It 
is then natural to ask whether 6 in the above statement can be improved to ||(5||^i/2. 

Theorem 1.8. There exist radial initial data uq for which ||'Uo||^i/2 < HQHi/i/a and 
u{t) blows-up in finite forward time. 



This follows from Theorems 1.4 and 1.5 by considering certain Gaussian initial 



data with negative phase (see Figure 6.3 ). It needs to be remarked, however, that this 
theorem relies on one piece of numerical information - that value ||Q||^i/2 = 27.72665. 
It is an analytical result in the sense that one need not numerically solve the NLS 
equation. 

This analytical result is further supported numerically for a variety of nonreal initial 
data with the inclusion of negative quadratic phase: Uq{x) = 0(x)e*''''^' with radial 
and real-valued and 7 < 0. We note, however, that we did not observe any real-valued 
initial data uq with ||Mo||iji/2 < ||Q||ifi/2 evolving toward finite-time blow-up solutions. 
Hence, we pose the following conjecture: 
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Conjecture 3. If the initial data Uq is real-valued and ||Mo||£fi/2 < ||Q||jifi/2, then u{t) 
scatters as t ^ — oo and t —>■ +00. 

Assuming Conjecture 3 holds, we elaborate further in Conjecture 4 below. When 
working with a one-parameter family of profiles, the H^^'^ norm can apparently predict 
both scattering and blow-up if the profiles are monotonic; this is summarized in our 
next conjecture. 

Conjecture 4. Consider a real-valued radial initial data profile ip{r) that is strictly 
decreasing as r +00. Let ao = ||<5||i/i/2/||^/'||£fi/2. Then the solution with initial 
data uo{x) = ai/j{x) scatters if a < ao and blows-up if a > oq. 

However, if the profile is not monotonic, then the H^^^ norm appears to only give a 
sufficient condition for scattering. This is illustrated in the simulations for oscillatory 



Gaussian data - see Figure 9.1 



1.1. NLS as a model in physics. 

1. Laser propagation in a Kerr medium [26l ^U\. This model is inherently two 
dimensional (in xy) and the time t in fact represents the ^-direction, and is derived 
via the paraxial approximation for the Helmholtz equation. The nonlinearity arises 
from the dependence of the index of refraction on the amplitude of the propagating 
wave. Blow-up in finite time is observed in the laboratory as a sharp focusing of 
the propagating wave. Ultimately, the non-backscattering assumption, and hence the 
NLS model, breaks down. 



2. Langmuir turbulence in a weakly magnetized plasma [331 ES]. A plasma is 
modeled as interpenetrating fluids of highly excited electrons and positive ions. The 
Langmuir waves propagate through the electron medium. The principle mathematical 
model is the Zakharov system |33], which is a nonlinearly coupled Schrodinger and 
wave system. The Schrodinger function is a slowly varying envelope for the electric 
potential and the wave function is the deviation of the ion density from its mean 
value. The NLS equation arises as the subsonic limit of the Zakharov system, which 
is obtained by sending the wave speed +00. Blow-up in finite time is the central 
phenomenon of study in |33], since it predicts the formation of a cavern of shrinking 
radius confining fast oscillating electrons whose collisions dissipate energy (at which 
point the model breaks down). The Zakharov model is inherently 3d, although certain 
experimental configurations can be modeled with the Id or 2d equations. 

3. Bose-Einstein condensate (BEC) [6J. BEC consists of ultracold (a few nK) dilute 
atomic gases where u gives the wave function (e.g. the number of atoms in a region E 
is and the coefficient of the nonlinear term is related to the scattering length 
hy g = Sna. The scattering length depends upon the interatomic potential and can 
be either positive or negative. While the model is inherently 3d, the imposition of a 
strong confining potential in one or two directions can effectively reduce the model to 
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two or one dimensions, respectively. Experiments showing blow-up are reported for 
®^Rb condensates in [7j and for ^Li condensates in [Ti] . 

In each situation, blow-up is physically observed (although of course the model 
breaks down at some point prior to the blow-up time). 

1.2. Acknowledgements. S.R. thanks Pavel Lushnikov for bringing to her attention 
his 1995 paper on the dynamic collapse criteria. J.H. and S.R. are grateful to Gadi 
Fibich for discussion and remarks on a preliminary version of this paper. S.R. is 
partially supported by NSF grant DMS-0808081. J.H. is partially supported by a 
Sloan fellowship and NSF grant DMS-0901582. 



2. Inequalities 



In this section, we prove the two inequalities (1.14) and (1.15) needed for Theorems 
OlandOl 



2.1. Inequality (1.14). Of course (1.14), the uncertainty principle, is standard, al- 
though we include a proof for completeness. By integration by parts, 

2 



^ y (V ■ x) c/x = Re j{ 



X ■ Vm) u dx . 



By Cauchy-Schwarz, 
9„ 

This provides the sharp constant since the inequality is achieved by only the above 
application of Cauchy-Schwarz, which is saturated when xu = CVu, i.e., when m is a 
Gaussian. 







2 




2 


+ 






/(x.V.)r.dx 





2.2. Inequality (1.15). In this section, we prove the following 



Proposition 2.1. The inequality 
(2.1) 



holds with sharp constant C 



1/14 



22.7='-7r^ 



1.3983. Moreover, all functions for 



which equality is achieved are of the form (3(j){ax), where 



\x 



2\l/2 



[X 



ifO<\x\<l 
if \x\ > 1. 



First we prove there exists some constant C for which (2.1) holds. Given u such 
that ||xm||2,2 < oo and < oo, define v by u{x) = av{Px) with a and jS chosen 

so that llxf 1^2 = 1 and \\v\\l4 = 1, i.e.. 



10 _ 3 

— 1 1 W 1 1 1 1 XW 1 1 7 



/3 



I W 1 1 ^4 1 1 ecu 1 1 
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By Holder, 



Ill2 ||||2 ||||2 

rllL2 = Ir IIl2(|x|<1) + Ir llL2(|a.|>i) 



4 \i/2 



a;|<l) + lk'^llL2(|2,|>i) 



3 



4 \i/2 



< 



Slj.j +1, 

which completes the proof of (2.1) with nonsharp constant C = ({^ttY^'^ + 1)^/^ ^ 
1.7455. 

Remark 2.2. We can optimize the above splitting argument by splitting at radius 
r = (t^j:) , which gives the constant C ~ 1.7265. However, this is still off from the 
sharp constant C ~ 1.3983 stated in Prop. 2.1 The reason is that the estimates 
applied after the splitting lead to an optimizing function which is the characteristic 
function of the ball of radius r. However, such a function fails to have both \\u\\l4, = 1 
and ||xti||L2 = 1. 



We now proceed to identify the sharp constant C in (2.1) and the family of opti- 
mizing functions. Consider the Lagrangian 



3 4 
I 7 II 7 



m 



defined on X =^ fl L'^{{x)'^dx), x G 



L2 



Lemma 2.3. 

(1) Any minimizer (j) (without loss of generality taken real-valued and nonnegative) 
of L in X is radial and (nonstrictly) decreasing. 

(2) There exists a minimizer cj) of L in X. 



Remark 2.4. Lemma 2.3 does not say that a minimizer (f) needs to be continuous, 
strictly decreasing, or compactly supported. We only discover this to be the case in 
the next step of the proof. 

Proof. We first argue that any G X can be replaced by a radial, monotonically 
(perhaps not strictly) decreasing function such that 

(2.2) ||01U. = ||0|U., 1101^4 = 1101^4, 



(2.3) ||x0|U2 < ||x0|U2 



Moreover, we have equality in (2.3) if and only if = 0, which occurs if and only if 
itself is radial, nonnegative, and (nonstrictly) decreasing. Given 0, let 

-E</,,A = {a; I |0(x)| > A } , ;U<^(A) = |E<^,a|- 
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Recall that (by Fubini's theorem) 

r+oo 

(2.4) 11011^,= / p\^-'ii^{\)d\. 



Jx=o 

Note that /i^ : [0, +cx)) [0, +oo) is a nonstrictly decreasing right-continuous func- 
tion. Thus, 

lini fi^{X) = /i^(A) < lim fi^a) 

(t\A (t/A 

with equality if and only if A is a point of continuity. Since fi^ is nonstrictly decreasing, 
there may be intervals on which /i^ is constant that we need to exclude in order to 
achieve invertibihty. Let be the set of all r such that Ai<^^({|7rr^}) has positive 
measure (is a nontrivial interval). Note that is an at most countable set (because 
fifP is nonstrictly decreasing). For all r ^ F, define the radial function on M.^ by 

4 

0(r) = A iff /i<^(A) < -vrr^ < lim fi^{a) . 

O (T/"A 



Note that n,f, = fx^, i.e., and are equidistributed. By (2.4), we have (2.2). Let 

h^{X) =W |a;|^ dx. 

By Fubini, 

r+oo 

11x011^2 =2 / Xh^{X)dX, 
Jx=o 

and similarly for 0. For each A, we have /i^(A) < /i<^(A), since = \E^^x\ but E^^ 

is uniformly positioned around the origin (it is a ball centered at 0). Consequently, 
(g holds. 

We note that the above argument establishes (1) in the theorem statement. To 
prove (2), we need to construct a minimizer by a limiting argument. Let 

m *== inf L((h) . 

Let 0„ be a minimizing sequence. By approximation and the above argument, we 
can assume that each 0„ is continuous, compactly supported, radial, nonnegative, 
and nonstrictly decreasing. By scaling (0„(x) ^ q;0„(/?x)) we can also assume that 
||0n||L2 = 1 and ||0n||L4 = 1- We have ||a;0n||L2 < ||x0i||i2 =^ A. Since for each n, the 
function (f)n{f) is decreasing in r, we have 

^'>INn|lW)= / ^^p'\Up)?dp>\Ur)?^. 

Jp=0 3 

i.e., |0n('")| ^ r~^/^. Similarly, working with the fact that ||0n||L4 = 1, we have that 
for all n, |0n('")| ^ r~'^/^ . Combining these two pointwise bounds, we have 

(2.5) |0„(^)|< ^-3/4(1 +^)-7/4^ 
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with implicit constant uniform in n. Thus, for each r > 0, the sequence of nonnega- 
tive numbers {(t>n{.f)}n=i is bounded and hence has a convergent subsequence. By a 
diagonal argument, we can pass to a subsequence of 0„ (still labeled (j)n) such that for 
each r G Q, we have that lim„^+oo 0n(^) exists. Denote by (j){r) the limiting function 
(for now only defined on Q). 

We claim that 0„ converges pointwise a.ej^ Pass to a subsequence (still labeled 0„) 
such that for each r G 2~"N, ^ < r < ra, we have \4>n{,f) — 0('")| < 2^^*^. Let 

C(0 = + 1)2"") - 2-2" for j2-" < r < (j + 1)2"", 

0+(r) 0(j2-") + 2-2" for j2-" < r < (j + 1)2-". 

Then for each n, 

0;(r) < 0„(r) < 0+(r) 

and 0- is a pointwise nonstrictly increasing sequence of functions, while 0+ is point- 
wise nonstrictly decreasing sequence of functions. Since is a decreasing function, 

(0(j2-") - 0((j + 1)2-")) = 0(l/n) - 0(n) < n''' 



On — loE n 



by (2.5). Thus, there can exists at most ?7,^/^2"/2 indices j for which the jump 
0(j2-") — 0((j + 1)2"") > 2-"/2. Thus, the measure of the set if„ on which 
0+(r) -0;(r) > 2-"/2 satisfies < n3/42"/22-" < This establishes 

that for a.e. r, 0(r) =^ lim^^+oo 0n('") exists and moreover for a.e. r, 

lim 0^(r) = 0(r) = lim 0^(r). 

n— >+oo n-^+oo 

We know that ||a;0-(x)|U2ci < ||a;0„|k2. By monotone convergence, we conclude 
that 

||x0(x)||i2 = lim ||x0;(x)||^2(i<^,<„) < lim \\x(pn\\L^ = m^l^ . 
Similarly, we have that 

|l4<1. 



By (2.5) and dominated convergence, we conclude that 



L2 = lim ||0„||l2 = 1 . 

n— >+oo 

Hence, 1^(0) < m (and thus, -^^(0) = m), i.e., G X is a minimizer. □ 



^Note that is a sequence of functions each of which is decreasing, but it is not the case that 
for each r, the sequence of numbers (/>n(r) is decreasing. Thus, proving that (/'n(r) converges for a.e. 
r > is a httle more subtle. 
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Lemma 2.5. If cf) E X , cf) is radial, nonnegative, (nonstrictly) decreasing, and solves 
L'{(f)) = 0, then there exist a>0, /3>0, 0<7<1 such that (j){x) = j3(j)^{ax), where 

for \x\ > 7. 



(2.6) 
Proof. Writing 



\x\ 



[X 



logL(0) = — log||x(/ 



li^ + ^iog 



2 

L2, 



it follows that for each ip, 

O = L'(0)(V^) = L(0) 



3|x| 



+ 



L2 

From this, we see that it must have the form 



4 



L2 



ip dx . 



(2.7) 



0(X) =/?(!- «2|^|2)l/2. 



on the set where (f){x) 7^ 0. Since (j){x) is decreasing, we see that (2.7) holds on 
< |a;| < 7/a for some < 7 < 1, and 0(a:) = for |x| > 7/a. □ 



Let G X be a minimizer for L{(j)) as in Lemma 2.3 We know that must solve 
the Euler-Lagrange equation L'{(j)) = 0, and hence, Lemma 2.5 is applicable, which 
establishes that (j){x) = P(j)^{ax) for some 7, < 7 < 1, where 0^ is given by (2.6). 
We now plug 0(x) = ^cp^^ax) into L{(j)) to determine 7. However, since L{(f)) is 
invariant under the rescaling 0(x) i— > P(j){ax), it suffices to take a = 1, /3 = 1 in this 
computation. We compute 



4 

L4 



5 7 



\X(, 



Atx 



7 



and we thus obtain 



m 



14 



2 

L2 



47r 



27" 



5 



(4vr) 



2 f 1 



A tedious computation shows that 7 = 1 produces the minimum value, which is 

35-52 



14 



22 . 75 . 7^2 



This completes the proof of Prop. 2.1 



3. New blow-up criteria 



In this section, we prove Theorem 1.4, 1.5, and 1.6 
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3.1. The blow up criteria of Lushnikov. Here, we prove Theorem 1.4 It is an 
adaptation and further investigation of the dynamic criterion for coUapse proposed 
by P. Lushnikov in [50]. 

We write M = M[u], E = E[u], etc., for simphcity. By the virial identity (1.4) in 
the case n = 3, 



(3.1) Vtt{t) = 24:E - 4:\\Vu 



and the bound (1.14), we obtain 

M2 l\Vtit)\^ 



V{t) 4 V{t) 

Now, rewritting the equation (3.2) to remove the last term with V^^ by making a 



(3.2) Vttit) <2AE-9 

Now, rewritting the equatic 
substitution V = B'^'^ , we get 

(3.3) B„<30£i3'/»-f|l. 

This differential inequality is an equality with some unknown non-negative quantity: 

(3.4) 5.<30E5^/^-|JJ-/(t). 



The equation (3.4) is the key in further analysis and in [20] is called the dynamic 
criterion for collapse. To analyze this equation, Lushnikov proposes to use a mechan- 
ical analogy of a particle moving in a field with a potential barrier. Let B = B{t) be 
the position of a particle (with mass 1) in motion under 2 forces: 

Btt = Fi + F2, 

where 

dU 225 
(3.5) Fi = -— with the potential U = -25EB^/^ + —M^B^^^, 

and 

F2 = —g'^it), some unknown force which pulls the particle towards zero. 

If this particle reaches the origin in a finite period of time, 

B{t*) = for some < t* < cx), 

then collapse necessarily occurs at some time t <t*. 
Several observations are due: 

• If the particle reaches the origin without the force —g'^{t) (i.e., if B{ti) = some 
< ti < 00), then it also reaches the origin in the situation when this force is 
applied (-8(^2) = for some < t2 < ti) . 

• If E < 0, then B always reaches the origin and collapse always happens. Thus, the 
more interesting case to consider is E > 0. 
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The potential U in (3.5) as a function of 5 is a convex function (at least for positive 
B) with the maximum attained at 

_/3M2y/^ _ , 75v^ M3 

Bmax — I o ~^=r I and Umax — U [B„ 



8 E J -K^raax, ^y^- 

Define the "energy" of the particle B: 
(3.6) s^t) = ^^ + Um)), 



which is time dependent due to the term g{ty in (3.4). However, recall that for our 
purposes it is sufficient for B to reach the origin if B satisfies only Bu = Fi (see the 
first observation above), for which the energy S{t) is conserved. 

The analysis is facilitated if we introduce the rescaled variables B{s) and S{s), 
where 

^, , ^ ~fl6Et\ ~fl6Et\ 16E 

Bit) = Bmax B j , Sit) = Umax£^-^) , ^ = ^t. 



We obtain 



If we set 



then (3.6) converts to 



16 V 53/5 



S{s) = -B,{sr + U{B{s)). 



The potential U{B) is depicted in Figure 3.1 From this energy diagram, we can 
identify two sufficient conditions under which B{s) necessarily reaches in finite time: 

(1) S(0) < 1 and B(0) < 1. In this case, the value of -8^(0) does not matter. 

(2) £^(0) > 1 and -8^(0) < 0. In this case, the value of -B(O) does not matter. 

Now define V = B^'^" . Then 



E = \v^'^{V^ - + 3) . 



Introduce the function 



(3.7) f{y)= jl- + v-^. 



V 



2 



We obtain 



£<l ^ \Vs\< f{V) 
S>1 ^ \Vs\>f{V) 
Thus, we see that sufficient condition (1) for blow-up above equates to 

V^(O) < 1 and - /(r(0)) < V;(0) < /(V^(0)) 
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2r 



1.5 



0.5 



-0.5 









1 








U(B) 







0.5 



1.5 



2.5 



Figure 3.1. A depiction of the rescaled potential U as a function of 
position B. In the case S < 1, a particle starting at 5(0) < 1 is trapped 
in that region and moves to the origin in finite time. In the case S > 1, 
a particle with initial velocity -Bs(O) < 1 has sufficient energy to reach 
the origin in finite time, regardless of its initial position B[0). 



and condition (2) equates to 



K(o) > -f{vm. 



This is graphed in Figure 3.2 



The two separate conditions can be merged into one: the solution blows-up in finite 
time if 

'+/(F(0)) if\/(0)<l 
-f{V{0)) ift>(0)>l 

Tracing back through the rescahngs, we see that the relationship with V(t) and V{s) 
is 



v;(o) < 



which completes the proof of Theorem 1.4 



3.2. An adaptation. In this subsection, we prove Theorem 1.5 The analysis here 



is similar to that used above in the proof of Theorem 1.4, except that we use (1.15) 



in place of (1.14). 
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Figure 3.2. The two sufficient conditions (1) and (2) described in the 



text ^3.1 for blow-up in finite time convert to the conditions on K(0) 



in terms of ^(0) depicted in this figure. 



By energy conservation, we can rewrite the virial identity as 

Vtt = 2AE - A\\Vu\\l, 

•8 

= 16E -2||m||^4. 



By (1.15) 



2M2 

(3.9) Vtt < IQE ^ 

Let 



7 



4M2 

Uiy) = -IQEV - 



Then, as in the previous subsection, we can interpret V{t) as giving the position of a 
particle subject to a conservative force —dvUiV) plus another unknown nonconser- 
vative force pulling V{t) toward 0. The corresponding mechanical energy is 

(3.10) s{t) = lv,' + U{V). 

Restricting to the case > 0, we compute that U {V) achieves its maximum Umax at 
Knax, with 

mI -12mIeI 
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To facilitate the rest of the analysis, we introduce a rescaling. Define V{s) and £{s) 
by the relations 

V(t) = VmaxV{at) , £{t) = \Umax\£{at) , a = 7 , s = at. 



Me 



Then 



Vssis) <l{l- V-'/'{ 



and (3.10) equates to 



2 ' 3 3 

From this, we identify two sufficient conditions for blow-up in finite time. 

(1) £{0) < -1 and V{0) < 1. 

(2) ^(0) > -1 and K(0) < 0. 

Let f{V) be defined by Q. Then 



S<-1 ^ \Vs\<\llfiV), 



£>-l ^ \Vs\>\lhiV). 



Thus, condition (1) above holds if and only if 

V{0)<1 and -^-f{vm<Vs{0)<^f{Vm 
and condition (2) holds if and only if 

K(o)<-y|/(no)). 

Clearly, we can merge the two conditions into one: we have blow-up in finite time 
provided 

V 3 |_/(V(0)) if V(0) > 1. 
Substituting back V{t), we obtain with to = 4:C^^/^E'^/^M~''^/''^V{0) 
Vt{0) C'l^ ^ f2 //H if < 1 



M 2V3(ME)V6 V3l-/(cj) ifo;>l. 



which gives (1.16) and finishes the proof of Theorem 1.5 
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3.3. Infinite variance radial case. Let ^/'(x) be a smooth, radial, (nonstrictly) 
increasing function such that 



|a;p if |x| < 1 
2 if Ixl > 2 



Let Vji denote the localized variance, defined in (1.17). Note that Vr/M < 2R^. We 



need to replace inequality ( |1.15[ ) with a localized version. 
Lemma 3.1. Suppose that Vr/M < . Then 



1 

2 2 TT \ 



4 3 



Proof. Let = 2Vr/M so that, by assumption, we have < B?. Then 

^ = Ikllia = ||M|||2(|^.|<.r) + ||m||^2(|^.|>^) 

^1/2 

< I - vrr^ ' 



/4 1 

4 „o M 



3 1 IL,I|2 



where we note that the inequality ||w|||2(|^|>r) < :^Vr requires < R^. The result 
follows. □ 

We have not made any effort to identify the sharp constant here. 



Now we turn to the proof of Theorem L6 By direct calculation, we have the local 
virial identity: 



Vliit) = 4 1 d.du^lj (I) d.udku ~ Ja^j (I) \u\' - / 
Note that 

V;^it) = WE - 2\\uit)\\i, + Aniuit)) , 

where 

AR{u{t)) < -^Ik||i2(|^|>^) + ||M||i4(|^|>K) . 

Recall that we are assuming ME > 1. Then, if we take 
(3.11) R^>S-^M\ 
we have 

1 11 S 

■^\m\Lmx\>R) <^<0E. 



\U 



2 
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Also, by the radial Gagliardo-Nirenberg inequality, 



6 \\'Vu\\l2 



Hence, 



< 

27T \\u\\l2 

< — ||Vm||l2E^/^ since ME > 1 
27r 

< \6E + \6\\Vufr2 
- 4 4 II 11^ 

< 6E + 6\\u\\l,. 

Vl^{t) < (16 + 26)E - (2 - 6)\\u\\l4 . 



Now we follow the analysis in §3.2 but we must ensure that for all t, 
(3.12) ^ < -R\ 



In fact, (3.12) will act as a bootstrap assumption that will be reinforced by the 
mechanical analysis. Suppose that (3.12) holds, and thus. Lemma 3.1 is applicable. 
Then 

VII < (16 + 25)E - (2 - 5)M'/'{C^r'V^'/\ where = (^^) ^ ■ 



Let 

UiVn) = -(16 + 25)EVr - 2(2 - 25)(Coo)"^M^/V^~'/' 



and define the mechanical energy as 

The maximum of UiVR) occurs at Vmax and is equal to Umax, where 

{1-6)ImI 6{8 + 6)1 {1- 6)1 mIeI 



14 
3 



+ 5)3(Coo)^S3 (a 

We introduce a rescaling: Define V^(s) and £{s) by the relations 

6i(8 + 5)l(Coo)5S^ 



V{t) = VmaxV{at) , £(t) = \Umax\ S{at) , « ... 

(1 - 6)3M 



1 , ,7 

6 



Then 



1 ~ n ~ ,~, ~ , ~ , Hnf 1 ~ 2~ 
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-0.9 r 



The maximum occurs at V = 1 



-1.1 



-1.2 ■ 



-1.3 - 



-1.4 



-1.5 - 



-1.6L 



0.5 



1.5 



V 



v = v, 

V 




Figure 3.3. A depiction of UiV), with H indicated. Note that 
V'(O) < Vb and always H > 1. 



The maximum of U{V) now occurs at \^ = 1. The bootstrap assumption (3.12) 
equates to 



(3.13) 



Vr{s) < 



R\8 + 6)-siC^)^E 



14 „2 
3 



2(1 - 5)3M3 

where we have defined the right-hand side as H, the "bootstrap threshold". But by 



(3.11) and the assumption ME > 1, we have 



Vb > 



(1-5)35 

We thus have Vf, 3> 1, provided 6 is taken sufficiently smallj^ 

The (rescaled) potential U{V), and Vb are depicted in Figure 3.3 From this, we 



identify two sufficient conditions for blow-up in finite time, noting that in each case. 



if (3.13) holds initially, then it will hold for all times: 



(1) ^(0) < -1 and V{0) < 1. 

(2) ^(0) > -1 and V;(0) < 0. 



The smallness on S here does not depend on M, E, etc. It depends only on ijj{x), the weight 
appearing in the local virial identity. 
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The remainder of the analysis is the same as in the previous section, and thus, The- 
orem 11.61 is estabhshed. 



Finally, we present the proof of Corollary 1.7 



Proof of Cor. 1.7 assuming Theorem 1.6. Take R = fiMS where /i is the con- 



stant in the second equation in (1.18). Since Uq is real, (1.19) converts to the state- 
ment 

(3.14) (M£)i'^# = M^«l 



(see the graph of g in Fig. 

Vr{0) 



1.2 



M3 
Decompose 



Ml 



1 



|x|<(5i/2M(M£;)-i/3 



iIj{x/ R)\uq\^ dx 



+ 



M 



|x|>5i/2Af(M£;)-i/3 



■iIj{x/R)\uo\'^ dx 



I + 11. 



For I, note that, when |x| < 5^''^ M{M E)-^/^ < R, 

ipix/R) = \x\yR^ = \x\'^6/iJ^M^ < 52(ME)-2/V^ 
and thus, |I| < S'^{ME)~'^^^jj,~'^. For II, just use \ip{x/R)\ < 1 and the assumption 



( [r20| to obtain |II| < 5^{ME)-^/^. Hence 

Vr{0) 



MR^ 



< 5\ME) 



-2/3 



Thus, the first condition in ( |1.18[ ) is satisfied, and moreover. 



M3 MR^5 



< f/S{ME) 



-2/3 



Therefore, (3.14) holds provided S is sufficiently small. 



□ 



4. PRELIMINARIES FOR THE PROFILE ANALYSES 



In the sections that follow (^|5j-[9j), we consider several different initial data families. 
Here we record some facts needed. 

4.1. H^^^ norm. For radially symmetric data Uq, the Fourier transform can be ex- 
pressed as 



poo 

uo{R) =2R^^ / uo(r) sm{2TTRr)rdr 
Jo 



(4.1) 



2nR-'/' 



Uo[r) Ji/ 



Ji/2{2nRr)r^^'^dr 
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where J 1/2 is the Bessel function of index | given by 



Ji/2(27ri?r) 



sin(2'7ri?r) 



We can then obtain 
(4.2) 



;7r2 / \uo{R)\^ dR. 
Jo 



4.2. Properties of Q. Recall that the Pohozhaev identities (1.6) hold, which in the 
case n = 3 take the form 

(4.3) l|VQ||i. = 3i|Q||i. and HQUi^ = 4||Q||i.. 

It then follows that E[Q] = \M[Q] 01 \\\VQ\\l. 



We computed numerically: 

IIQIli2(]R3) = 



wvQr 



iiQiiiv. 



An I Q^{r)r^dr ^ 18.94, 
Jo 

POO 

An / g2(r)r^c/r ^ 20.32, 
Jo 

\Q{R)\^R^ di? 27.72665 . 



4.3. Reformulation of blow up conditions. For real valued initial data. Theorem 

simplifi 

ViO)< 



L4| condition (1.12) can be simplified to 

3M2 
8^^ 



(4.4) 



Similarly, the condition (|1.16|) from Theorem 1.5 can be simplified to 

1 



(4.5) 



V{0) < c with c 



4 (714/3 



and C from (1.16). 



Observe that the second condition (4.4) is an improvement over (4.5) for real valued 
initial data when 



5^2 



M[u]E[u] > ^M[Q]E[Q] ^ 2.06 M[Q]E[Q]. 
45U 



Thus, when discussing the real valued initial data, we will refer to (4.4) and (4.5) 



instead of Theorems L4 and L5, correspondingly; the blow up conditions look simple 
in this case. 



For complex valued initial data. Theorem 1.4 condition (1.12) has to be considered 



separately for positive and negative values of \4(0). Define (a scale invariant quantity) 

8EV{0) 



3 M2 

then for positive valued Vt{0) blow up happens 



BLOW-UP CRITERIA FOR 3D CUBIC NLS 29 

• in the intersection of the regions: 

(4.6) 0<.<1 and ^ < 2^ f ^J^'l^'^.^^^sf, 

and for the negative valued Vt{0) blow up happens 

• in all of the region < uj < 1 

• and in the intersection of the regions: 

In this paper we consider the complex valued initial data with the quadratic phase 

(4.8) «o(r) = /(r)e^^^^ 

with /(r) - real valued radial function and 7 G M. For such initial condition we have 

POO 

V{0)=AttF, Vt{0) = 32TT-fF with F= r^|/(r)|2dr. 

Jo 

Also 

(poo poo \ 

2tt J r^lV/l^c/r-vry r'^\f\'^dr\+8TT-f'^F = E^ + E\ 

Since 

(4.9) Vt{0)^ - 32V{0)E'-' = 0, 



the second condition in (4.6) is simplified further to 

[3 M 8E^V(0) 

and analogously, the second condition in (4.7) will be as above inequality with the 
reversed sign (we write it out for future reference) 

Similarly, Theorem 1.5 condition (1.16) has to be studied separately for positive 
and negative values of V^(0). Denoting (also a scale invariant quantity) 

we obtain that for the positive valued Vt{0) blow up happens 

• in the intersection of the regions: 
(4.12) 

, ^^(0) 2V2 f 4Ci4/W(0) 
'<^^' V<C^^[cVWm^ ^^ -3(ME)V3 

and for the negative valued Vt(0) blow up happens 



1/2 
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• in all of the region < k < 1 

• and in the intersection of the regions: 
(4.13) 

|V^(0)| . 2^2 f M'^l^ 



K> 1 and 



M 



> 



For the initial data with the quadratic phase as in (4.8), the second condition in (4.12) 
reduces to 



(4.14) 



M3/2 



C7F(0)i/2 



+ 4 



EV(0) fME\ 



1/3 



M2 



-3 



> 0, 



due to (4.9). The second inequality in (4.13) reduces to the same inequality as above 



except with the reversed sign (again, we write it out for the convenience of future 
reference): 



(4.15) 



h4 ^-3 <0 



In computations below we study the conditions (4.6)-(4.7), (4.10 )-(4.11 ) and (4.12) 



(4.15) instead of (1.12) and (1.16), respectively. In graphical presentation we refer 



to the set of conditions (4.6)-(4.7), (4.10)-(4.11 ) as "Condition from Thm. I.4" and 



(4.12) - (4.15) as "Condition from Thm. 1.5 



5. Q PROFILE 

In this section we study initial data of the form 
(5.1) Uo{r) = \"^Q{\r)e'^'\ 



where Q is the ground state defined by (1.5). Note that Uq has been scaled so that 
M[u] = M[Q] for all A > 0. 



We compute 



E[Q] 



3A2 - 2X^ + 



37 
A2 



||VMo|li2 

WW 



.2 , 7 



A^ + 



L2 



A2 



where, by a Pohozhaev identity (|4.3|), 
(5.2) 



4 J 



, WyQWh 



livgf 



4 JyQWh 

o7 



L2 



WQf 



L2 



The blow up when energy is negative occurs when 

2 



(5.3) 



f < 



A^ 
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< 1 and 



> 



L2 



For positive energy blow-up is analytically proven in the case 

E[Q\ 

1 by Theorem 1.1 These conditions equate to 

A'(l - A') < f < - SA' + 2A=^) = ^(1 - A)'(l + 2A) , 

which is not possible for any < A < 1. On the other hand, scattering is proved in 
the case ^^j^ < 1 and pjQ\\2'^ < 1; which just reduces to the condition < A < 1 and 

f <-A2(1-3A2 + 2A3). 
3 



This is depicted in Figure 5.1 



0.25 



0.2 



0.15 



0.1 - 



0.05 




2A3) 

l|v^jl^><4yQ||2 



0.2 0.4 0.6 0.8 1 

A 

Figure 5.1. The curves for the mass-gradient and mass-energy thresh- 



olds for the Q profile in (5.1) 



We now investigate the conditions for blow up from Theorems |1.4| and |1.5[ Note 
that 

v{Q) = ^A\yQ\\l and vm = i^A\yQ\\l- 



A 



A2 
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Hence, depending on the sign of 7, Theorems |1.4[ |1.5| provide different ranges of 



values (A, 7) for which blow up occurs. We start with 7 < 0. In what follows we use 



7 instead of 7, recall their simple relation (5.2) 



Theorem 1.4 we investigate the region given by < < 1, see notation in (4.6) 



and then its complement with an additional restriction (4.11). The condition u < 1 
is 



(5.4) 



f < + 



1 WQW 



M\yQ\\l 



note that the right side is nonnegative when A > 



(4.11) is 



4||yQiii 



) ^ 1.15. The condition 



(5.5) 



f > -A^ 
' 3 



WQWl 



WyQWl 



4(iA 



\\yQ\\] 



+ 3 



- 1 



provided A > 



A^ 



3 iigii 



4||yQII 



.45. 



Thus, for 7 < 0, or negative Vt{0), union of (5.4) and (5.5) will give a region where 



solution will blow up in finite time. It turns out that the first condition (5.4) is a 



part of the region covered by (5.5), and the last one is depicted in Figure 5.2 under 



the name "Blow up by Thru. I.4 



For 7 > 0, the blow up region by Theorem 1.4 is an intersection of (5.4) and the 



inequality (5.5) with the reversed sign, which is depicted in Figure 5.3 under the name 



"Blow up by Thru. \l-4\ '- It turns out that this region has previously been covered by 
our Theorem ll.il 

Summarizing, Theorem L4 provides a nontrivial blow up condition for the pair 
(A, 7^) only if 0.45 < A < 1.15 and 7 < (see Figure Is^. 



Theorem |1.5[ similarly to the above, we start with 7 < and investigate the 



region given by < k < 1, see notation in (4.12), and then its complement with an 



additional restriction (4.15). The condition k < 1 is 



(5.6) 



f < 



5-3^/^ ligil^ 

S7r-75/2 ll^/Qli; 



+ 



A^-A^ 



with the right side being nonnegative when A > 



condition (4.15) is 

(5.7) f>tx'-\' + -i^X^ 



3iiQir 



1 WQWl 
C'hQh 



1.25. The 



\\yQ\\l\ . , JyQW 



WQWl 



A + 2- 



2\ 3 
2 \ 



WQWl 
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where C is from (1.16), and is valid for any A > 0. Again, it turns out that the first 



condition (5.6) is a part of the region covered by (5.7), and the last one is depicted 



in Figure 5.2 under the name "Blow up by Thru. 1.5". 



For 7 > 0, the blow up region by Theorem 1.5 is an intersection of (|5.6) and the 

under the 



inequality (5.7) but with the reversed sign, which is depicted in Figure 5.3 
name "Blow up by Thm 



1.5". This region has also been previously covered by our 



Theorem II. 11 

In summary. Theorem L5 provides a nontrivial blow up condition for the pair 
(A, 7^) for any < A < 1.25 and 7 < (see Figure 5.2). In comparison with Theorem 
|1.4[ it provides a wider range for < A < 0.762. 



Note that Theorems 1.4 and 1.5 provide new information on the blow up behavior 



for the Q profile of the form (5.1) with negative phase. In particular. Theorem 



1.4 



gives a nontrivial range of 7 when 0.45 < A < 1.15. Further extension is given by The- 
orem 1.5 for any A < 0.762, see Figure 5.2 Both Theorems provide a blow up range 



"under" the mass-gradient condition, thus, showing that the last condition is irrele- 
vant for determining long time behavior in the region when M[-u]i?[M] < M[(5]-E'[Q]. 



Lastly, we compute ||mo||^i/2 norm using (4.2) numerically and then compare all 
conditions about the global behavior for this initial data together with numerical 



data in Figures 5.2 and 5.3 with negative and positive signs in the initial phase. 



correspondingly. For the positive phase Theorems L4 and L5 do not provide any new 
information, however, we have numerical range on blow up threshold and include the 
plot for illustration and completeness. 



Example 5.1. Consider uq = Q{r)e'^'^^ (i.e., take X = 1 in (5.1)). Note that 

WyQWr 



M[uo]E[uo] 



1 + 



WQWl 



M[Q]E[Q] > M[Q]E[Q], 



and thus, Theorems Ll_ and 1.3\ can not be applied. However, the solution with such 
initial condition will blow up in finite time if 

\ 1/2 



7 < 



IIQII^ 



V 



WyQWli^- 



4 iiQii; 



3 WyQW] 



WQWl 
M\yQWl 



-0.177, 



by Theorem 1.4, or when 



7 < 



WQWi 



54C^\\yQ\\l \\\yQ\\ 



WQW2 ^2C'\\yQ\\ 



2\ 3 
2 



WQW 



WQWl 
M\yQWl 



1/2 



-0.279, 



by Theorem L5. In this example. Theorem \l.4\ is more powerful than Theorem L5 



however, this is not always the case, as can be seen from Figure 5.2 (for example, for 



A < .762 Theorem 1.5 gives a larger range for blow up). 
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Figure 5.2. Global behavior of the solutions to the Q profile initial 
data with the negative quadratic phase: Uo{r) = A^/^(5(Ar) e"*'"^'''^. 
Here, 7 is the renormalized 7, see (5.2), namely, 7^ ~ 1.437^. The 
region of "theoretical scattering" is provided by Thm. see also 

Figure 5.1 "Blowup by Thm. 1.4' is given by (5.5) and "Blowup by 



Thm. 1.5' is given by (5.7). The intersection of these two conditions 
occurs at A ~ 0.762. 



5.1. Conclusions. 



is not valid; the nu 
curve in Figure 



5.2 



(1) The condition "||tio||£fi/2 < ||(5||^i/2 implies scattering" 
merical blow-up curve is below the ||'Uo||^i/2 = IIQII^i/2 
This supports Conjecture 2. 

(2) The condition "HmqIIl^ || V-UqIIl^ < IIQIU^ || V(5||l2 implies scattering'^ is not 



valid (unless M < M[(5]-E'[Q] as in Theorem |1.1 ); not only the nu- 
merical blow-up curve is below the ||mo||l2|| Vmo||l2 = ||(5||l2 || VQHl^ curve in 



Figure 5.2, but also both Theorems 1.4 and 1.5 provide range of (A, 7) for 



which blow up from the initial data 5A_ occurs and this range is below the 
||mo||l2||Vmo||l2 = ||Q||l2||VQ||l2 curve in Figure 1^2, 



(3) Previously, no theoretical blow-up result for the profile (5.1) with < A < 1 



could be obtained from Theorem 1.1 The new blow-up criteria in Theorems 
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Figure 5.3. Global behavior of the solutions to the Q profile initial 



r e 



-j|7|r^ 



. As 
The 



data with the positive quadratic phase: Ui^ir) 
before, 7 is the renormalized 7, see (5.2), namely, 7^ ~ 1.437 
region of "theoretical scattering" is the same as in Fig. 5.2 and is pro- 
vided by Thm. 1J_ "Blowup by Thm. |1.4f is given by the complement 
of (5.5) intersected with (5.4) and "Blowup by Thm. 1.5' is given by 
the complement of (5.7) intersected with (5.6). 



1.4 and 1.5 give a nonempty set of (A, 7) with 7 < for which blow up occurs. 



see conditions (5.4)-(5.5) and (5.6)-(5.7) as well as the illustration in Figure 

E2 



6. Gaussian profile 
In this section, we study initial data Mq of the form 

(6.1) Mo(x) = pe-"^''/2 e'^"' , r = |x| , x G Ml 
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By scaling, it suffices to consider the cases 7 = (real data) and 7 = The main 
parameters are 

„2 37r3/2p2 / 2 

l|Vno||i.-^(l + 42, 
7^ \ 



M[u] = 




a3/2 ' 


E[u] = 


7r3/2p2 / 


V{0) = 


3 7r3/2p2 


2a5/2 ' 


II l|2 27rp2 



,7' 



1277r3V 



a5/2 

2\ 1/2 

1 + 4-' 



a;2 



To compute the last expression ||no||^i/2, consider the Fourier transform of Uq (here, 

i?2 = |^|2^ ^ e R') 



uq{R) = p I — e ^^+47^ e ^^+47^ , 



a — 



2^77 

where we used e^(''^'+2^^) dz = A/^e"^^"/", a, 6 G C. By Q we have 



6.1. Real Gaussian. Take 7 = in (6.1). Then by scaling, the behavior of solutions 
is a function of p/ ^/oi. We have 

• E[u] > if 

(6.3) p < \fa ^ 2.91 \fa\ 

• the condition on the mass and gradient Hwollia || V'Uo|||2 < ||Q||l2 || V(5||i2 implies 

(6.4) p < 2^/^ 71-'/' IIQIU2 ^ ^ 2.19 Va; 

• the mass-energy condition M[-u]i?[M] < M[(5]-E[(5] is 



which gives 

(6.5) p < 1.92^/a and p > 2m^/a. 
• the invariant norm condition ||uo||^i/2 < ||Q||^i/2 is 

(6.6) p < (27r)-i/2 ligil^,/, ^ ^ 2.10v^. 
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• (Theorem 1.4) the condition (4.4) is 

(6.7) p > (iVtj 

• (Theorem 1.5) the condition (4.5) is 



1/2 



a ^ 2.38^/a 



{6.i 



3 . 23/2 . 75/2 X 1/2 
',307ri/2 + 75/2 ) 



Numerical simulations: the results for the real Gaussian initial data (6.1) are in 



Table |6.1[ For p > pt the blow up was observed, for p < Ps the solution dispersed 
over time. For example, for a = 1 the threshold is between 2.07 and 2.08. This 
is consistent with the previously reported threshold by Vlasov et al. in |29j {p = 
2.0764). From this table it also follows that ps/y/a G (2.07,2.075) and Pb/y/ci G 
(2.077,2.08). 



a 


0.5 


1.0 


2.0 


4.0 


6.0 


8.0 


10.0 


Ps 


1.46 


2.07 


2.93 


4.15 


5.08 


5.87 


6.56 


Pb 


1.47 


2.08 


2.94 


4.16 


5.09 


5.88 


6.57 



Table 6.1. Thresholds for blow up/scattering from numerical simula- 
tions for the Gaussian initial data: for p < Ps scattering was observed 
and for p > pb blow up in finite time was observed. For comparison the 



values of p from (6.6) are also listed. 



In Kuznetsov et. al. [23] it is reported that for the Gaussian initial data (6.1) 
with 7 



(6.9) 



and 



(6.10) 



0, the numerical condition for collapse is when two conditions hold: 

|Vmo||1I|mo||2 



llvgiiiiiQii; 



> 0.80255 or 



(2 ■0.80255)^/4.,^., 

P> ._3/4 WQWl' 



a 



2.0759 Wa 



E[uo]M[uo] 
E[Q]M[Q] 



> 1.1855 or 2.0764V« <p< 2.6105Va. 



Thus, the numerical threshold p/ 
consistent with this report, see also Figure [6T 



2.0764 is reported in [23]. Our data is 



To compare the conditions (6.4) - (6.8) with the numerics, we graph them in Figure 



6.1 For clarity of presentation, and also for comparison with the 7 7^ case considered 
— p 

next, we plot —= on the vertical axis vs a on the horizonal. 

a 
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E = 



I 



V(0) = cM'''E-2/3 



V(0) = %M'E-i 



l|VUoll2l|Uoll2=l|VQ|l2l|Q|l2 



E<0 



M(u)E(u)<M(Q)E(Q) 
& 

l|VUoll2l|Uoll2>l|VQ|l2l|Q|l 



> V(0) < %M^E-' 




M(u)E(u) = M(Q)E(Q) 



numerical threshold for 
blowup & scattering 



M(u)E(u) < M(Q)E(Q) 

r & 

l|VUoll2l|Uoll2<l|VQ|l2l|Q|l2 



Figure 6.1. Global behavior of the solutions with the (real) Gaussian 
initial data u^ir) = pe"^'/^^ The line denoted by 1/(0) = S/SM^E"^ 
is the threshold for blow up from Theorem 1.4, see (6.7); similarly, the 
line denoted by V^(0) = cM^/'^i?"^/^ is the threshold for blow up from 
Theorem 1.5, see (6.8). The line "theoretical scattering" is given by 
Thm. |1.1[ see (6.5). The numerical threshold (dashed line) comes from 



Table 



6.1 normalized by ^/a. For all other values refer to text in ^6.1 
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6.2. Gaussian with a quadratic phase. Now we consider (6.1) with 7 7^ 0. (By 
scahng it suffices to consider 7 = ±|.) We compute 

• E\u] > if 



(6.11) 



P 



< (6 V2) + 4X^ ^ 2.91 (1 + 4^ 



2\ 1/2 



the condition on the mass and gradient H^^ollia || Vuo||^2 < ||Q||l2 II V(5||i2 is 



(6.12) 



p < 



2V4||g|| 



L2 



a 



2.19 Va 1 + 4 



7r3/4 (1 + 4^)1/4 
the mass-energy condition M[M]ii^[M] < M[Q]i?[Q] is 



,7 



a' 



-1/4 



(6.13) 



4q;2 



3 1 + 4 



a" 



a 



1 

2' 



< i^llQllt- 



The left hand side is a cubic polynomial in y"^ = p'^/a and can be solved explicitly 
as a function of a, though with a very complicated expression. We 



to obtain 



p 



1/2 



list a few values for 7 = ±1 in Table 6.2 The inequality in (6.13) holds for p < pi 



a 


0.25 


0.5 


1 


2 


4 


6 


8 


10 


Pi 


0.41 


0.79 


1.45 


2.42 


3.70 


4.62 


5.38 


6.04 


P2 


6.01 


4.60 


4.09 


4.46 


5.65 


6.75 


7.71 


8.58 



Table 6.2. The positive real roots of the equation in (6.13). 



and p > p2. 

the invariant norm condition ||uo||^i/2 < ||Q||^i/2 is 

1/2 



(6.14) 



p < a 



27.72665 



2 7r(a2 + 472)1/2 



2.10 Va 1 + 4 



a' 



-1/4 



(Theorem 1.4) the condition uo <1 from (4.6) amounts to 



(6.15) 



P 



> (4^2)^^^ (^1 + 6 



7 



a' 



2 \ 1/2 



2.38A/a 1 + 6 



,7 



2\ 1/2 



a" 



similarly, uj >1 from (4.7) will be the above with the reversed sign. The condition 



(4.10) for positive 7 with y = pj \foL is 

/ - 6v^ + 4^^ + 64 > 



(6.16) 



and with the reversed inequality sign for negative 7. The last inequality is cubic 
in y^ producing two positive roots, which are listed in Table 6.3 for 7 = ±|. 
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a 


0.25 


0.5 


1 


2 


4 


6 


8 


10 


Pb 


0.45 


0.86 


1.58 


2.68 


4.22 


5.37 


6.33 


7.16 


Pt 


6.01 


4.60 


4.08 


4.39 


5.42 


6.35 


7.17 


7.91 



Table 6.3. The positive real roots, pt and pb, of the function in (6.15) 



for the condition (4.10) (Theorem 1.4). 



(Theorem 1.5) the condition k < 1 from (4.12) amounts to 

1/2 



2\ 1/2 
2 



similarly, k > 1 from (4.13) is the above inequality with the reversed sign. The 



condition (4.14) with y = pl\f(x for positive 7 is 

1/3 



(6.18) 



15V27r 



V 



1/3 



> 0, 



75/2 23/2/^ \2-7^ J V\ «V 23/2 

and with the reversed sign for negative 7. The positive real zeros of the function 



in (6.18) for 7 = ±^ are listed in Table 



6.4 



a 


0.25 


0.5 


1 


2 


4 


6 


8 


10 


Pb 


0.48 


0.93 


1.68 


2.81 


4.39 


5.57 


6.55 


7.41 


Pt 


6.01 


4.61 


4.10 


4.46 


5.54 


6.51 


7.36 


8.13 



Table 6.4. The positive real roots, pt and pb, of the function in (6.18) 



for the condition (4.14) (Theorem 1.5) 



Numerical simulations: we fix the quadratic phase 7 = ±.5 and vary the parameter 
a in our numerical simulations in order to obtain the blow up threshold, see data 



in Table 6.2 Note that we obtain different thresholds depending on the sign of 7. 



To compare conditions (6.11) - (6.18) with the numerics, we graph the dependence 
of p on a in Figures j6.2 and separately for positive and negative values of the phase 
'-' — "^^^^ on the vertical axis to observe the asymptotics as a — > 00, and 



7 



±i. We plot 



thus, approaches the case of the real Gaussian initial data. 

6.3. Conclusions. The above computations show 

(1) Consistency with Conjecture 3: if Uq is real, then ||mo||^i/2 < ||(5||jji/2 implies 
u{t) scatters (numerical data). 

(2) Consistency with Conjecture 4 (since uq is based upon a radial profile that is 
monotonically decreasing): if ||Mo||£fi/2 > ||Q||ii-i/2, then u{t) blows-up in finite 
time (numerical data) . 
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a 


0.5 


1.0 


2.0 


4.0 


6.0 


8.0 


10.0 


pt 


2.8 


3.00 


3.56 


4.58 


5.43 


6.17 


6.83 


pt 


2.8 


3.01 


3.57 


4.59 


5.44 


6.18 


6.835 


Ps 


0.79 


1.44 


2.42 


3.76 


4.76 


5.59 


6.313 


Pb 


0.80 


1.45 


2.43 


3.77 


4.77 


5.60 


6.315 



Table 6.5. Thresholds from numerical simulations for the Gaussian 
initial data with phase 7 = ±.5: for p > pt blow up in finite time 
was observed and for p < pf scattering was observed; + superscript 
indicates the threshold for the positive phase 7 = .5 and — superscript 
indicates the negative phase 7 = —.5. 



(3) Theoretical proof of the Conjecture 2, i.e., "||mo||^i/2 < ||Q||^i/2 implies scat- 
tering" is false for an arbitrary radial data. We show that it is possible to pro- 
duce a radial Gaussian initial data with negative phase (e.g., for < A < 3.3 
by Theorem 1.4) such that ||Mo||£fi/2 < IIQII^fi/2 and the solution u{t) blows 
up in finite time. 

(4) The condition "||mo||l2 || Vmo||l2 < ||Q||l2 II^QIU^ implies scatterincf^ is not 
valid (unless M[u\E[u] < M[Q]E[Q] as in Theorem 

(5) In all three cases (real data, data with positive phase and with negative phase) 
Theorems lA and L5 provide new range on blow up than was previously 
known from our Theorem II. II 

(6) For the Gaussian initial data with negative phase and small values of a (0 < 
a < 2, see Figure 6.3) the numerical threshold for scattering and blow up 
coincides with the scattering threshold provided by Theorem 1.1 As a — 00 
the numerical threshold (for both positive and negative phases) approaches 
the values ||mo||^i/2 = ||Q||j^i/2. 



7. Super-Gaussian profile 
Now we modify the initial Gaussian data to the "super Gaussian" profile: 

(7.1) Mo(r) =pe-"'^*/2e'^"', r = \x\, x eR\ 
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S_E^O M(u)E(u) < M(Q)E(Q) 

y & 

l|VU„||2||U„||2>||VQ||2||Q||2 



Condition from Thm. 1.5 
Condition from Thm. 1.4 



IVl{u)E(u) = IV1{Q)E(Q) 

numerical threshold for 
blowup & scattering 



M(u)E(u)<M(Q)E(Q) 



VUoll2l|Uo||2<||VQ||2||Q||2 



Figure 6.2. Global behavior of the solutions to the Gaussian initial 
data with the positive quadratic phase uo{r) = pe~°''^ /2e+*2'' . The 
curve denoted by "Condition from Thm. 1.4| ' comes from the largest 



positive root of the equation in (6.15), namely, values Vtl \foi in Table 



6.3 Similarly, the curve denoted by "Condition from Thm. |1.5[' comes 



from the largest positive root of the equation in (|6.18|), namely, values 
Ptl \/ot in Table 



6.4 



The curve "theoretical scattering" is provided by 
Thm. see (6.13) and values pi in Table 6.2 The numerical thresh- 
old (dotted curve) comes from values p^, p'l in Table 6.2 normalized 

by 



'a. 
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Figure 6.3. Global behavior of the solutions to the Gaussian initial 
data with the negative quadratic phase uo{r) = pe~"^' /^e~*2'' . The 
curve denoted by "Condition from Thm. |1.4[' comes from the smallest 



positive root of the equation in (6.15), namely, values Pb/\/c( in Table 



6.3 Similarly, the curve denoted by "Condition from Thm. |1.5[' comes 



from the smallest positive root of the equation in (6.18), namely, values 
Pb/y/ci in Table 6.4 The curve "theoretical scattering" is the same as 
in Fig. 6.2 and is provided by Thm. |1.1[ see (6.13) and values pi in 
Table 6.2 The numerical threshold (dotted curve) comes from values 
p~ , p^ in Table 6.2 normalized by y/a. This plot illustrates that both 



Theorems L4 and |1.5| provide ranges of the Gaussian initial data Uq 
such that ||Mo||jji/2 < ||(5||^i/2 and u(t) blows up in finite time. 
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For this initial data we calculate 

.2 /g^ 



M\u\ 



E\u] 



Here, Vis + 11 



a3/4 V 4 



5 + 4 



7 



4y2aV4r(|) 
2y2a5/4r(|) 



TT 5 + 4 



7 



[r (1)]^ P 

21/4 
2y27rV7 



7.1. Real super Gaussian. When 7 = 0, we have 
• > if 



(7.2) 



p < 



r(!) 



(7.3) 



the condition on the mass and gradient ||mo||l2|| VMo||i2 < ||<51|l2||^QIIl2 implies 

23/8 3I/4 



p < 



\\Q\\L^a'/'^2.10a'/'; 



51/47^3/4 

the mass-energy condition M[M]ii^[M] < M[Q]ii^[Q] is 



4^2, 



57r 



[r(!)] 



3M2 p2 



a 



21/4 



<2llQlli^' 



and thus, we obtain 



(7.4) 



P 



<1.71a^/^ and j9>3.44a^/^ 



(7.5) 



(Theorem 1.4) the condition (4.4) is 

21/8^1/2 



p > 



5 - 



6 



1 4' 



1/2 



a 



(7.6) 



p > 



a 



(Theorem 1.5) the condition (4.5) is 

2I/8 . 5I/2 . 7^3/4^7/2 

r(|)(w/2 + 4[r(|)]4)i/2 

the invariant norm condition ||mo||^i/2 < HQ II ^1/2 is given in Table 
first compute the Fourier transform of Uq using (4.1) and then H^^^ norm by (4.2) 
which is listed in the second row of the Table 17.11 for various a. The third row 
indicates the values of p, denoted by pi/2, for the threshold in the invariant norm 
condition. We observe that Pi/2/v^ ~ 2.02. 



7.1 



Here, we 
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Numerical simulations: the results with the super Gaussian initial data (7.1) are in 
Table 7.1 For p > Pb the blow up was observed, for p < Ps the solution dispersed 
over time. We observe that ps/ ^fa ^ 2.01 and Pb\^ ~ 2.02. For convenience Table 



7.1 contains the if^/^ norm calculations. 



a 


.25 


.5 


1 


2 


4 


6 


8 


10 


^1 Uo\\hi/2 


13.54 


9.58 


6.77 


4.79 


3.39 


2.76 


2.39 


2.14 


Pl/2 


1.43 


1.70 


2.02 


2.41 


2.86 


3.17 


3.41 


3.60 


Ps 


1.42 


1.69 


2.01 


2.39 


2.85 


3.15 


3.39 


3.58 


Pb 


1.43 


1.70 


2.02 


2.40 


2.86 


3.16 


3.40 


3.59 



Table 7.1. The H^^^ norm of the super Gaussian initial data depend- 
ing on a and values of p for the H^^'^ condition as well as the numerical 
results for blow up threshold and global existence - p > Pb and p < Ps, 
correspondingly. 



7.1 



To compare the conditions (7.2) - (7.6) with numerical data, we graph the depen- 
dence of p on a in Figure 



For clarity of presentation we plot —= on the vertical 



P 



axis. 



7.2. Super Gaussian with a quadratic phase. When 7 7^ 0, we have 
• E[u] > if 

91/8 1/2 / 2\ 1/2 / ^2 " 

p<^^(5 + 4^) ai/^^1.58ai/^ (5 + 4^ 
r(|) V «/ \ a 

the condition on the mass and gradient ||%||i2 1| VuoHia < ||Q||l2 II V(5||i2 implies 



(7.7) 



03/80I/2 / 2\ -1/4 

(7.8) p < ||g|U2 a^/^ ^ 3.15 a^/^ 5 + 4— 



the mass-energy threshold M[M]i?[M] = M[Q]ii^[Q] is 



(7.9) 



4v^i 



a 



TT 5 + 4 



7 



a 



p 



21/4 



a 



7.2 



the real positive zeros of the above expression when 7 = ±| are in Table 
for the invariant norm condition ||mo||^i/2 < IIQII^i/2) we compute the Fourier 



Transform of Uo{r) = pe "'"■*/2 gn''^ again by (4.1) and then H^^^ norm by (4.2). 



The values of p when the condition ||Mo||i/i/2 = ||<5||iji/2 are in Table 7.3 as well. 
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a 


.25 


.5 


1 


2 


4 


6 


8 


10 






Pi 


1.00 


1.28 


1.60 


1.96 


2.37 


2.64 


2.85 


3.02 








3.34 


3.48 


3.81 


4.32 


5.01 


5.49 


5.88 


6.20 




Table 7.2 


. The positive 


real roots of the equation in 


7.9). 


a 


.25 


.5 


1 


2 


4 


6 


8 


10 


^1 Mo|Lyl/2 


17.53 


11.03 


7.29 


4.97 


3.45 


2.80 


2.42 


2.16 


Pl/2 


1.26 


1.59 


1.95 


2.36 


2.83 


3.15 


3.39 


3.59 



Table 7.3. The H^^^ norm of the super Gaussian initial data with 
the phase 7 = ih| depending on a is listed in the second row and values 
of p for this condition are in the last. 



• (Theorem 1.4) the condition u < 1 from (4.6) is 
(7.10) 

4\ 1/2 



P > 



21/8^1/2 



5 + 4 



7 



a 



6 



r 



a 



1/4 



1.58 3.63 + 4 



7 



2\ 1/2 



a 



a 



and the condition ( 4.10[ ) (with y = p/a^^^) is 
57r2 



(7.11 



TT 



6[r(f)]^ 



25/43[r(|)]2 ^ + 



23/23l/2[r(|)]2 



vr 5 + 



Mill! ,,2 
21/4 y 



17^ >0' 



the positive real zeros of the left hand side with 7 = ±^ are in Table 



7.4 



a 


.25 


.5 


1 


2 


4 


6 


8 


10 


Pb 


1.62 


2.04 


2.56 


3.18 


3.90 


4.38 


4.75 


5.05 


Pt 


3.32 


3.43 


3.71 


4.13 


4.69 


5.10 


5.42 


5.68 



Table 7.4. The positive real roots, pt and p^, of the function in (7.11 ) 
for the condition (4.10) (Theorem 1.4). 



(Theorem 1.5) the condition (4.12) is 

21/8^3/4 (J7/2 



P > 



(7.12) 



r(|) (cv/2 + 4[r(|)]4) 

2\ 1/2 



1/2 



1+4 



r 



2 \ 1/2 



a 



a 



1/4 



1.29 1 + 4 



r 

a 



a 



1/4 
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and the condition ( |4.14[ ) with y = p/a^^^ is 
(7.13) 
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25/4[r(|)] 



3^12 / y + 



3 ■ 2^6 



5 + 4^ 
a I 



'3m2 



m)] 6 

-r 



1/3 



> 0, 



the positive real zeros of the left hand side with 7 = ±^ are in Table 



7.5 



a 


.25 


.5 


1 


2 


4 


6 


8 


10 


Pb 


1.41 


1.84 


2.36 


2.50 


3.68 


4.15 


4.51 


4.80 


Pt 


3.30 


3.41 


3.66 


4.05 


4.59 


4.97 


5.27 


5.53 



Table 7.5. The positive real roots, pt and pb, of the function in (7.13) 



for the condition (4.14) (Theorem 1.5). 



Numerical simulations: the results for the super Gaussian with the quadratic phase 

The rows denoted by p^, p^ 



7 



are given in Table 



7.6 



and p^ are 



thresholds for blow up and global existence/ scattering for the positive phase 7 = +7 



and negative phase 7 



2' 



correspondingly. 



To compare the conditions (7.7) - (7.13) with the numerics, we graph the dependence 

in Figures 7. 
phase 7 = ±|. We plot 



of p on a in Figures 7.3 and 7.2 separately for positive and negative values of the 

p — 

—= on the vertical axis to observe the asymptotics as a — 00 



a 


.25 


.5 


1 


2 


4 


6 


8 


10 


pt 


2.18 


2.27 


2.47 


2.76 


3.14 


3.41 


3.63 


3.81 


pt 


2.19 


2.28 


2.48 


2.77 


3.15 


3.42 


3.64 


3.82 


PJ 


1.02 


1.33 


1.69 


2.12 


2.60 


2.92 


3.17 


3.38 


Pb 


1.03 


1.34 


1.70 


2.11 


2.61 


2.93 


3.18 


3.39 



Table 7.6. Numerical results for blow up threshold {p > pt) and 
global existence {p < pf). The sign in the superscript indicates the 
positive or negative sign of 7, correspondingly. 



Remark 7.1. Numerical simulations showed that, for example, when a = 2 for 
2.3995 < p < Q the blow up occurs over the origin. When p > 10 the blow up 
happens on a contracting sphere. This phenomena we originally discussed in our 
heuristic analysis in [T^, it was also numerically obtained in [llj . 
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E = 



V(0) = %M^E-i 



V(0) = cM^/3E-2/3 



L E<0 I M(u)E(u)<W 
J1IVU0II2IIU0II2 



M(Q)E(Q) 

2> IIVQII2IIQII2 



> V(0) < cM'/3E-2/3 



IIVUoll2l|Uoll2=l|VQ|l2l|Q|l: 



M(u)E(u) = M(Q)E(Q) 



numerical threshold for 
blowup & scattering 



i 



IVI(u)E(u)<M(Q)E(Q) 
& 

||VUoll2llUoll2<l|VQ|l2l|Q|l2 



Figure 7.1. Global behavior of the solutions with the super Gaussian 



initial data UQir) = pe' 



-ar'^/2 



The line denoted by V{0) = S/SW^E' 



is the threshold for blow up from Theorem 1.4, see (7.5); similarly, the 
line denoted by ^(0) = cM'^/^E~'^/^ is the threshold for blow up from 
Theorem 1.5, see (7.6). The line "theoretical scattering" is given by 



Thm. 

Table 



1.1[ see (7.4). The numerical threshold (dashed line) comes from 
7.1 normalized by ^fa. For all other values refer to text. 
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Condition from Thm. 1.4 
Condition from Thm. 1.5 




IVI(u)E(u) < IV1(Q)E{Q) 
> & 

VUoll2l|Uoll2> IIVQII2IIQII2 



numerical threshold for 
blowup & scattering 



M(u)E(u)<M(Q)E(Q) 
& 

VUo||2||Uo||2<||VQl|2l|Q||2 



Figure 7.2. Global behavior of solutions to (1.1) with the super 
Gaussian initial data with the positive quadratic phase Uo{r 



pe 



The curve denoted by "Condition from Thm. 



1.4 



comes from the largest positive root of the equation in (7.11), namely, 



values pt/ \/a in Table 7.4 Similarly, the curve denoted by "Condition 
from Thm. |1.5|' comes from the largest positive root of the equation in 



(7.13), namely, values ptj -sfoL in Table 7.5, The dotted line "numerical 
threshold" is plotted from the values / ^fa. and pi / ^/a (indistinguish- 



able) from Table 7.6 The curve "theoretical scattering" is provided by 



Thm. |1.1| see (7.9) and values pi in Table 7.2 
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Figure 7.3. Global behavior of solutions to (1.1) with the super 
Gaussian initial data with the negative quadratic phase Uo(r) = 



The curve denoted by "Condition from Thm. 



1.4 



pe 

comes from the smallest positive root of the equation in (7.11), namely, 
values pb/ in Table 7.4, Similarly, the curve denoted by "Condition 
from Thm. |1.5[ ' comes from the smallest positive root of the equa- 
tion in (7.13), namely, values Pb/\^ in Table 7.5 The dotted line 



"numerical threshold" is plotted from the values p~l j \foL and p~^ j ^fot 
(indistinguishable) from Table 7.6 The curve "theoretical scattering" 
is the same as in Fig. 7.3 and is provided by Thm. |1.1| see (7.9) and 



values p\ in Table 7.2 
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7.3. Conclusions. The above computations show 

(1) Consistency with Conjecture 3: if uq is real, then ||uo||^i/2 < ||(5||jji/2 imphes 
u{t) scatters. 

(2) Consistency with Conjecture 4 (since Uq is based upon a radial profile that is 
monotonically decreasing): if ||-Uo||£fi/2 > IIQII^fi/2, then u{t) blows-up in finite 
time. 

(3) The condition "||mo||l2 || VwoULa < ||(5||l2 || VQHl^ implies scattering" is not 



valid (unless M[u]E[u] < M[Q]E[Q] as in Theorem [Llj). 
(4) In all three cases (real data, data with positive phase and with negative phase) 
Theorems L4 and |1.5 provide new range on blow up than was previously 



known from our Theorem LI Furthermore, Theorem L5 provides the best 
results. 



8. Off-centered Gaussian profile 
Next we consider the off-centered Gaussian profile: 
(8.1) uo(x) =pr2e-"''%^^"', r=\x\,xeR^. 

For this initial data we calculate 



M[u] 

E[u]- 

V{0) 
II ii 



" 32^a7/2 

37r3/V 
64\/2a5/2 

105 TT^/ V 

12872^9/2^ 



|VMo||r2 — ^ 

' 32v^a5/2 



11 + 35 



r 

a' 



11 + 35 



315p2 
2i%¥ 



a- 

105 7 TT^/ V 
1672^9/2 



327r' 



R 



dR 



3V (1 + 2^) 
(1 + S)V2- 



8.1. Real off-centered Gaussian. When 7 = 0, we have 
• Elu] > if 



(8.2) 



p < 



32 f llV2^ 
~35 



1/2 



a 



3/2^ 7.11 a3/2. 



the condition on the mass and gradient ||Mo|li2|| V'Uo||^2 < ||Q||i2 1| V(5||i2 implies 



(8.3) 



p < 



||g|h«3/^^3.46a^/2; 



1657r'^ 

the mass-energy condition M[M]ii^[M] < M is 

457r^p*^ 



212^6 



ii-^)<-IMI 



4 

L2, 
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and thus, we obtain 
(8.4) p<2.73a^/2 and p>7Ma^/^. 

I TmrQT'iQ'ni' Ti/^Tm r'/^n ri i i r^n 11-7/^1, 

1^1/2 



the invariant norm condition ||uo||^i/2 < ||Q||^i/2 is 



.5) 



p < 



(3 7r)i/ 



(Theorems 1.4) the condition (4.4) is 

32 



(8.6) 



p > 



21 



/ 62^2 Y 



a 



^ 6.38 



(Theorems 1.5) the condition (4.5) is 
17) 



25.21/4.72.111/2 



p > 



3. 51/2(75 + 25. 33/2. 5%l/2)l/2 



a 



^ 5.93 a^/l 



Numerical simulations: the results for the off-centered Gaussian initial data (8.1) 
are in Table 8J^ For p > pb the blow up was observed, for p < Ps the solution 
dispersed over time. We obtain that psjo?^'^ ~ 3.57 and phjo?^"^ ~ 3.58. 



a 


.25 


.5 


1 


1.5 


2 


3 


4 


5 


Ps 


0.44 


1.26 


3.57 


6.57 


10.11 


18.5 


28.4 


40.0 


Pb 


0.45 


1.27 


3.58 


6.58 


10.12 


18.6 


28.5 


40.1 



Table 8.1. The numerical results for blowing up {p > pb) and global 
existence {p < ps) for the off-centered Gaussian initial data depending 
on the parameter a. 



To compare all the above conditions (8.2) - (8.7) with the numerics, we graph the 



dependence of p on a in Figure 
vertical axis. 



6.1 



P 

For clarity of presentation we plot on the 



8.2. Off-centered Gaussian with quadratic phase. When 7 7^ 0, we have 
• E[u] > if 



p < 



32 / V2 

35" 



1/2 



11 + 35 



2\ 1/2 



a 



3/2 ^ 2.14 (11 + 35 



2\ 1/2 



a 



3/2. 



the condition on the mass and gradient H^tollia || Vuo||/,2 < HQHia || VQ||^2 implies 



.9) p < 



11 \ 1/4 



157r' 



IIQII2 11 + 35 



a' 



-1/4 



a 



^ 6.29 (11 + 35 



)— 

a' 



-1/4 



a 



3/2. 



BLOW-UP CRITERIA FOR 3D CUBIC NLS 



l=- 



theoretical 
scattering 

I 



E = 



}| M(u)E(u)<M(C 
E<0 L & 
J l|VUoll2l|Uoll2> 



M(Q)E(Q) 

IIVQII2IIQII2 



- V(0) = %M2E-i 



V(0) = cM^/3E-2/3 



V(0) < cM'I^E 



M(u)E(u) = M(Q)E(Q) 



numerical threshold for 
blowup & scattering 



l|VU„||2||Uoll2=l|VQ||2l|Q|I 



i 



J_ 



M(u)E(u)<IVl(Q)E(Q) 
& 

l|VUoll2l|Uoll2< IIVQII2IIQII2 



Figure 8.1. Global behavior of the solutions to (1.1) with the real 



off-centered Gaussian initial data (8.1) 
3/8M^£'~^ is the threshold for blow up from Theorem 
similarly, the line denoted by 1^(0) = cM'^^'^E''^/^ is the threshold for 



The line denoted by V{fi) 
see 



1.4 



8.6) 



blow up from Theorem |1.5[ see (8.7). The line "theoretical scatter- 
ing" is given by Thm. |1.1| see (8.4). Observe that the numerical 



threshold (dashed line, values are given in Table 8.1) is away from the 
line 1 1 uo 1 1 ^1/2 = ||(5||£fi/2, see ( |8.5[ ), which is different compared to real 
Gaussian and real super Gaussian initial data, see Fig. QA_ and [7]T 
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the mass-energy condition M[M]ii^[M] < M[Q]ii^[Q] is 



.10) 



3^4 



457rV 



11 + 35^- ^ 



a" 



2iO\/2, 



IIQIl2<o, 



the real positive zeros of the left hand side when 7 = ±^ given in Table 



8.2 



a. 


.25 


.5 


1 


1.5 


2 


3 


4 


5 


Pi 


0.17 


0.65 


2.3 


4.58 


7.31 


13.84 


21.54 


30.27 


P2 


3.29 


5.14 


9.51 


15.14 


21.9 


38.26 


57.8 


80.05 



Table 8.2. Real positive zeros of the function in (8.10) when 7 = ±^ 



the invariant norm condition ||^^o||^i/2 < HQ II ^1/2 is 



(8.111 



^ ^ a^/^ ^ 3A3- 4^ 



+ 234)1/2 



;i+ 25)^2 



a 



3/2 



(Theorems 1.4) the condition (4.4) is 
32 ■ 2^/4 



(8.12) 



p > 



^ + ^V'^^/^ - 12.68 fo.253+^' 



and the condition (4.10) with y = p/a^^'^ is 

2^^/105 3 • 72 



(8.13) 



7a/21i 11 + 



357^ 



315^21/2 



32 



15^2 



>0, 



the positive real roots of which when 7 = ±1 are in Table 



8.3 



a 


.25 


.5 


1 


1.5 


2 


3 


4 


5 


Pb 


0.60 


1.79 


5.37 


10.25 


16.18 


30.63 


47.97 


67.80 


Pt 


3.29 


5.14 


9.49 


14.99 


21.50 


37.00 


55.27 


75.91 



a 



3/2 



Table 8.3. The positive real roots of the function in (8.13) when 7 



±^ for the off-centered Gaussian initial data. 



(Theorems 1.5) the condition (4.5) is 



p > 



25 . 21/4 . 72 



(8.14) 



3. 51/2(75 + 25. 33/2. 5%V2)i/2 



11 + 35 



11 



2\ 1/2 



a 



3/2 



1.79 11 + 35 



11 
a' 



2\ 1/2 



a 



3/2 
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^.15) 



(26 ■ 33/2 ■ 5V/2 - 7^) 
25/2 . 



-V 



2^ ■ 73 ■ 11 
33-52 



-847r |^2"(11 + 35^) - 'il^^Ply^ Z/^ > , 
the positive real zeros of which with 7 = ±| is in Table 8 



a 


.25 


.5 


1 


1.5 


2 


3 


4 


5 


Pb 


0.30 


1.16 


4.15 


8.46 


13.79 


26.94 


42.92 


61.04 


Pt 


3.29 


5.14 


9.45 


14.80 


21.05 


35.75 


52.97 


72.38 



Table 8.4. The positive real roots of the function in (8.15) when 7 



±^ for the off-centered Gaussian initial data. 



Numerical simulations: the results for the off-centered Gaussian initial data with 



nonzero phase (8.1) are in Table 8.5 For p > pb the blow up was observed, for 



p <Ps the solution dispersed over time. 



a 


.25 


.5 


1 


1.5 


2 


3 


4 


5 


pt 


n/a 


n/a 


5.78 


8.93 


12.68 


21.5 


31.9 


43.6 


pt 


n/a 


n/a 


5.79 


8.94 


12.69 


21.6 


32.0 


43.7 


Ps 


.18 


.69 


2.50 


5.11 


8.337 


16.2 


25.8 


36.8 


Pt 


.19 


.70 


2.51 


5.12 


8.338 


16.3 


25.9 


36.9 



Table 8.5. The numerical results for blow up and scattering thresh- 
olds for the off-centered Gaussian initial data with the phase 7 = ±^: 
blow up if p > pt and global existence/scattering ii p < pf. 



To compare all the above conditions (8.8) - (8.15) with the numerical data, we 



graph the dependence of p on a in Figures 8.2 and 8.3 For clarity of presentation we 

p — — 

plot -—J- on the vertical axis. 

8.3. Conclusions. The above computations show 

(1) Consistency with Conjecture 3: if Uq is real, then ||mo||^i/2 < ||(5||jji/2 implies 
u{t) scatters. 

(2) Consistency with Conjecture 4: Numerical simulations show that the blow up 
threshold line is higher than the line ||Mo||ifi/2 = ||Q||iji/2, although this profile 
is not monotonic. 
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Figure 8.2. Global behavior of the solutions for the off-centered 
Gaussian initial data with positive quadratic phase. The curve de- 
noted by "Condition from Thm. 1.4| ' comes from the larges t po sitive 



root of the equation in (8.13 ), namely, values pt/ in Table 
ilarly, the curve denoted by "Condition from Thm. 1.5| 



8.3 



Sim- 
comes from the 

largest positive root of the equation in (8.15), namely, values pt/\/oi 
The dotted line "numerical threshold" is plotte d fro m the 

8^ The 



in Table 8.4 



values / ^fa^ and / ^fa' (indistinguishable) from Table 
curve "theoretical scattering" is provided by Thm. |1.1| see (8.10) and 
Table lO 
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Figure 8.3. Global behavior of the solutions for the off-centered 
Gaussian initial data with negative quadratic phase. The curve de- 
noted by "Condition from Thm. |1.4| ' comes from the smalles t po sitive 
root of the equation in (8.13), namely, values ph/ \fa in Table 8.3 



Sim- 
comes from the 

■2 



ilarly, the curve denoted by "Condition from Thm. 1.5| 



smallest positive root of the equation in (8.15), namely, values p^l \Jot 
in Table 8^ The dotted line "numerical threshold" is plotte d fro m the 
values p~ / ^fa^ and pl j ^fa (indistinguishable) from Table 8.5, The 
curve "theoretical scattering" is provided by Thm. |1.1| see (8.10) and 
Table lO 
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(3) The condition "||mo||l2 || Vmo||l2 < ||(5||l2 || VQlUa implies scattering" is not 
valid (unless M[M]£^[n] < M[(5]ii^[Q] as in Theorem 1.1), for example, when 
there is a negative initial phase present. 

(4) In all three cases (real data, data with positive phase and with negative phase) 
Theorems L4 and |1.5 provide new range on blow up than was previously 

Furthermore, Theorem |1.5| provides the best 



known from our Theorem |1.1 
results. 



9. Oscillatory Gaussian profile 



Lastly we consider a Gaussian profile with oscillations (referred from now on as 
'oscillatory Gaussian'), it is a sign changing profile: 



(9.1) 



u^fi{x) =p cos(/?r) e ^ e*''"' . 



Here, we fix the Gaussian e itself and change the frequency of oscillation /5. These 
data generate a solution to NLS denoted by Ujs. We obtain: 

3/2 2 3/2 2 



AV2 



4^/2 

3/2 2 

^3/2 2 3/2 2 

V{0) = —^v{(3), Vt{0)=^—^v{(3), 



I6V2 



U/3,o\\'iri/2 = 327r^j9^ / R 



2V2 

) 

2 

r cos(/5r) sin(27ri?r) e"'' dr 



dR 



TX 



■p- 



where 



m(/?) = l + (l-/?2)e-/^V2^ 
.(/?, 7) = 3(1 + 7^) (3(1 + 7^) + 67^) + /5V) e-^^/^ 

6(/3) = (3 + (1 - W')e-'' + 2(2 - P')e-^'/') , 

^;(/?) = 3 + (3-6/?' + /?^)e-^'/^ 

2 r ^2 

erf(x) = I e ^ dt. 



We list some values of the H^^^ norm of upfi in Table 



9.1 
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/5 





0.25 


0.5 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


1\Wp,0\\%i/2 


7T 


3.046 


2.788 


2.101 


1.879 


2.901 


3.933 


4.922 


5.906 


Pl/2 


2.97 


3.02 


3.15 


3.63 


3.84 


3.09 


2.65 


2.37 


2.17 


P 


7.0 


8.0 


9.0 


10.0 


15.0 


20.0 


25.0 




l\Wf3,o\\%i/2 


6.890 


7.875 


8.859 


9.844 


14.765 


19.687 


24.609 


Pl/2 


2.01 


1.88 


1.78 


1.68 


1.37 


1.19 


1.06 



Table 9.1. The H^^'^ norm of U/j^o and values of p for which 

||M/3,o||^1/2 = IIQII^l/2- 

9.1. Real oscillatory Gaussian. First, we consider the real oscillatory Gaussian 
data (7 = 0). We have 

• E[u] > if 



(9.2) 



p < 



m 



the condition on the mass and gradient 1 1 M/3 1 1 1,2 ||V'U/3^o Ilia < IIQIlL^II^QIlia implies 



(9.3) 



p<nQh 



6 



1/4 



^7r3m(/3)a(/3,0), 
• the mass-energy condition M[M]i?[M] < M gives 

(9.4) 7r3m(/5) (a(/3,0) -/6(/3))p4_32||g||4 < 0, or p<pi and p > pa, 

where pi and p2, the real positive zeros of the cubic polynomial (in p^) above, given 
in Table 








0.25 


0.5 


1.0 


1.45 


2.0 


2.5 


2.77 


3.0 


3.61 


4.0 


Pi 


2.72 


2.76 


2.88 


3.27 


3.48 


3.16 


2.72 


2.55 


2.43 


2.22 


2.12 


P2 


3.81 


3.87 


4.06 


5.00 


6.84 


10.49 


13.11 


13.42 


13.31 


12.97 


13.14 



/3 


5.0 


6.0 


7.0 


8.0 


9.0 


10.0 


15.0 


20.0 


25.0 


Pi 


1.91 


1.76 


1.64 


1.53 


1.45 


1.38 


1.13 


0.98 


0.88 


P2 


14.75 


17.17 


19.81 


22.48 


25.17 


27.87 


41.47 


55.13 


68.82 



Table 9.2. The values of p in the mass-energy threshold for the real 
oscillatory Gaussian. 

the values of p, denoted by pi/2, for which ||m/3,o||^i/2 = HQ 11^^1/2 5 given in Table 
I9T1 
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• (Theorem 1.4) the condition (4.4) is 
(9.5) 



(9.6) 



(Theorem 1.5) the condition (4.5) is 



a(/5,0)l/2^;(/3)3/4 



(7r3/2m(/3)7/2 + 8C7 y(/5)3/2 6(/3))i/2 ' 
Numerical simulations: the results for the (real) oscillatory Gaussian initial data 



are in Table 9.3 For p > Pb the blow up was observed, for p < Ps the solution 
dispersed over time. 








0.25 


0.5 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


Ps 


2.93 


2.98 


3.11 


3.601 


4.583 


4.13 


3.29 


2.92 


2.72 


Pb 


2.94 


2.99 


3.12 


3.605 


4.585 


4.14 


3.30 


2.93 


2.73 




7.0 


8.0 


9.0 


10.0 


15.0 


20.0 


25.0 




Ps 


2.61 


2.53 


2.48 


2.455 


2.42 


2.49 


2.52 


Pb 


2.62 


2.54 


2.49 


2.456 


2.43 


2.50 


2.53 



Table 9.3. Numerical simulations for the (real) oscillatory Gaussian. 
Here, the blow up was observed for p > ph and global existence for 

P < Ps- 



9.1 



To compare all conditi ons (|9.2[) - (|9.6|) with numerical data, we graph the depen- 
dence of p on /3 in Figure 



We plot p/a/1 + on the vertical axis. 

9.2. Oscillatory Gaussian with quadratic phase. When 7 7^ we have 
• E[u] > if 

'a(A7)^'/' 



(9.7) 



p < 



m 



the condition on the mass and gradient ||tt/3||i2 1| Vu/3^o||i2 < ||Q|lz,2|| VQ||^2 implies 



(9. 



P<2iig||2 



6 



1/4 



7r3m(/3)a(/3,7), 
• the mass-energy condition M[M]i?[M] < M[Q]ii^[Q] gives 

(9.9) 7r=^m(/3) (a(/3,7) -p26(/3))p^-32||g||^ < or p < p} and p > pi, 
where p\ and pi, the real positive zeros of the polynomial above with 7 = ±|, 



given in Table 9.4 



BLOW-UP CRITERIA FOR 3D CUBIC NLS 



61 



/9 





0.25 


0.5 


1.0 


1.6 


2.0 


2.4 


2.66 


3.0 


4.0 


PJ 


2.42 


2.47 


2.62 


3.13 


3.39 


3.09 


2.74 


2.56 


2.39 


2.10 


7 
P2 


4.46 


4.50 


4.61 


5.29 


7.97 


10.94 


13.37 


13.92 


13.77 


13.40 



f3 


5.0 


6.0 


7.0 


8.0 


9.0 


10.0 


15.0 


20.0 


25.0 


7 
Pi 


1.90 


1.75 


1.63 


1.53 


1.45 


1.38 


1.13 


0.98 


0.88 


7 
P2 


14.95 


17.34 


19.95 


22.61 


25.28 


27.97 


41.54 


55.18 


68.86 



Table 9.4. The values of p for the mass-energy threshold for the os- 
cillatory Gaussian with the phase 7 



the values of p, denoted by py^^ for which ||«/3,o|| 



1^1/2 



IIQII %i/2^ are given in Table 



f3 





0.25 


0.5 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


^ '"0,/3|^i/2 


2 


3.384 


3.042 


2.165 


1.933 


2.959 


3.942 


4.922 


5.906 


7 

Pl/2 


2.81 


2.86 


3.02 


3.58 


3.79 


3.06 


2.65 


2.37 


2.17 




7.0 


8.0 


9.0 


10.0 


15.0 


20.0 


25.0 




^lh^0,/3 1^1/2 


6.890 


7.875 


8.859 


9.844 


14.765 


19.687 


24.609 


7 

Pl/2 


2.01 


1.88 


1.78 


1.68 


1.37 


1.19 


1.06 



Table 9.5. The values of ||m/3,o||^i/2/p^ for the oscillatory Gaussian 
initial data with phase 7 = ±| and values of py^ ^01 which the if^/^ 
norm threshold holds. 



• (Theorem 1.4) the condition (4.4) is 

(9.10) p>b{[3r'/'(a{f3,^)- 



2\ -1/2 



and the condition (4.10) is 
(9.11) 



{ai(3,0)-p'biP))-3) {a{(3,^)-p%if3)f\ipy/' > 0, 



the real positive zeros of which are in Table 9.6 



• (Theorem 1.5) the condition (4.5) is 
(9.12) p > 



23/4^7/2^3/4^(^^^)1/2 



7r3/2m(/3)l/2 + 23/2C7t;(/3)3/25(/5))l/2' 
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/3 





0.25 


0.5 


1 


1.6 


2 


2.4 


2.66 


3 


4 


P 


2.68 


2.74 


2.93 


3.83 


7.02 


9.86 


12.05 


12.58 


12.53 


12.55 


P 


4.39 


4.42 


4.52 


5.07 


7.88 


10.89 


13.30 


13.83 


13.67 


13.30 






5 


6 


7 


8 


9 


10 


15 


20 


25 


P 


14.29 


16.79 


19.47 


22.19 


24.91 


27.64 


41.31 


55.01 


68.73 


P 


14.86 


17.26 


19.88 


22.55 


25.23 


27.93 


41.51 


55.16 


68.85 



Table 9.6. Zeros of the function in the inequahty (9.11) Theorem 1.4 



and the condition (4.14) is 
(9.13) 



2i/2C7y(/?)i/2 m{py )^ m(/?)2 
277r3p4 



8C 



14 



m(/5) (a(/5,7)-/&(/5)) > 0, 



the positive real zeros of which are hsted in Table 9.7 








0.25 


0.5 


1 


1.6 


2 


2.4 


2.66 


3 


4 


Pb 


2.81 


2.87 


3.06 


3.90 


6.71 


9.10 


10.45 


10.57 


10.35 


10.56 


Pt 


4.46 


4.49 


4.60 


5.14 


7.91 


10.86 


12.05 


13.37 


12.01 


12.52 




/3 


5 


6 


7 


8 


9 


10 


15 


20 


25 


Pb 


12.22 


14.41 


16.74 


19.11 


21.49 


23.88 


35.92 


48.03 


60.15 


Pt 


13.96 


16.11 


18.44 


20.80 


23.18 


25.58 


37.61 


49.72 


61.85 



Table 9.7. The positive real zeros of the polynomial in (9.13), by 
Theorem 11.51 



Numerical simulations: the results for the oscillatory Gaussian initial data with the 

For p > Pb the blow up was observed, for p < Ps 



±1 are in Table 



9.8 



phase 7 

the solution dispersed over time. 



To compare all the above conditions (9.7) - (9.13) with the numerical data, we 



graph the dependence of p on a in Figures 9.3 and 9.2 For clarity of presentation we 

p — — 

plot — , on the vertical axis. 

9.3. Conclusions. The above computations show 

(1) Consistency with Conjecture 3: if uq is real, then ||Mo||/fi/2 < ||Q||hi/2 implies 
u(t) scatters. 
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/5 





0.25 


0.5 


1.0 


2.0 


3.0 


4.0 


5.0 


pt 


3.56 


3.60 


3.70 


4.07 


4.94 


4.75 


3.70 


3.19 


pt 


3.57 


3.61 


3.71 


4.08 


4.941 


4.76 


3.71 


3.20 


Ps 


2.42 


2.47 


2.63 


3.20 


4.263 


3.64 


2.98 


2.70 


Pb 


2.43 


2.48 


2.64 


3.21 


4.26 


3.65 


2.99 


2.71 


p 


6.0 


7.0 


8.0 


9.0 


10.0 


15.0 


20.0 


25.0 


pt 


2.92 


2.76 


2.66 


2.59 


2.55 


2.48 


2.53 


2.56 


pt 


2.93 


2.77 


2.67 


2.60 


2.56 


2.49 


2.54 


2.57 


PJ 


2.56 


2.47 


2.42 


2.38 


2.377 


2.37 


2.44 


2.48 


Pb 


2.57 


2.48 


2.43 


2.39 


2.378 


2.38 


2.45 


2.49 



Table 9.8. Threshold for blow up and scattering in numerical simu- 
lations for the oscillatory Gaussian with the phase 7 = ±|. 

(2) Consistency with Conjecture 4: observe that the real oscillatory Gaussian 
initial data uq is a radial profile that is not monotonically decreasing, thus, 
the condition ||mo||^i/2 > 1^1/2 does not necessarily imply that u(t) blows- 
up in finite time. 

(3) The condition "llwolUa || VmoIIl^ < IIQIU^ || V(5||l2 implies scatterinff^ is not 



valid (unless M[M]i?[M] < M[(5]-E'[(5] as in Theorem 1.1) even for the real 
oscillatory Gaussian initial data. 
(4) In all three cases (real data, data with positive phase and with negative phase) 



Theorems 1.4 and 1.5 provide new range on blow up than was previously 



known from our Theorem |1.1[ For small oscillations (small P) Theorem 1.4 



provides the best range for blow up and for large oscillations Theorem 1.5 
provides a better result. 
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M(u)E(u)<M(Q)E(Q) 

l|VUoll2l|Uoll2> IIVQII: 



M(u)E(u) = M(Q)E(Q) 



M(u)E(u)<M(Q)E(Q) 

J ||VU„||2||U„||2< IIVQII2IIQII2 



Figure 9.1. Global behavior of the solutions to (1.1) with the real 
oscillatory Gaussian initial data 9.1, with the rescaled vertical axi s 
p/ The blow up threshold curve denoted "by Thm. 



1.4 



is given by (9.5) and the blow up threshold curve denoted "by Thm. 
1.5' is given by (9.6). Observe that these curves intersect at ~ 1.23. 



Therefore, for small oscillations {(3 < 1.23) Theorem 1.4 provides the 
best range for blow up, correspondingly, for large oscillations (/3 > 1.23) 
Theorem 1 1 . 5| gives a better range. The curve "theoretical scattering" is 
given by Thm. |1.1[ see (9.4) and values pi in Table 9.2 Observe that 
the numerical threshold (dotted curve, values are given in Table 9.3) 
for small oscillations coincide with the curve ||Mo||/fi/2 = ||(5||£fi/2, see 
Table 9.1, and for large oscillations these curves separate. 
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M(u)E(u)<M(Q)E(Q) 
& 

l|VUoll2l|Uoll2> IIVQII; 



M(u)E(u) = M(Q)E(Q) 



M(u)E(u)<M(Q)E(Q) 

J l|VUoll2l|U„||2< IIVQII2IIQII2 



Figure 9.2. Global behavior of the solutions to (1.1) with the oscil- 
latory Gaussian initial data with positive quadratic phase 7 = +|, see 



(9.1). The vertical axis is rescaled p/yl + to compare with the real 
case. The blow up threshold curve denoted "by Thm. 



1.4 



by (9.11), see values pt in Table 9.6 



noted "by Thm. 1.5' is given by (9.13), see values pt in Table 9.7 



IS given 

the blow up threshold curve de- 



Observe that these curves intersect at (3 
oscillations (/3 < 1.89) Theorem 



1.4 



1.89. Therefore, for small 
provides the best range for blow 



up, correspondingly, for large oscillations {(3 > 1.89) Theorem 1.5 gives 
a better range. The curve "theoretical scattering" is given by Thm. 
|1.1| see (9.9) and values pj. The numerical threshold (dotted curve) 
is given in Table 9.8, values pf and p'^. All p values in this graph are 



normalized by ^yT+]p 
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M(u)E(u)<M(Q)E(Q) 

l|VUoll2l|Uoll2> IIVQII: 



M(u)E(u) = M(Q)E(Q) 



M(u)E(u)<M(Q)E(Q) 

J l|VUoll2l|U„||2< IIVQII2IIQII2 



Figure 9.3. Global behavior of the solutions to (1.1) with the oscil 

1 



latory Gaussian initial data with negative quadratic phase 7 



2' 



see 



(9.1). The vertical axis is rescaled p/yl + to compare with the real 
The blow up threshold curve denoted "by Thm. 



case. 



1.4 



IS given 



by the complement of (9.11), see values pb in Table 9.6 the blow up 



threshold curve denoted "by Thm. 1.5' is given by the complement of 



(|9.13|), see values pb in Table 9.7[ Observe that these curves intersect at 

pro- 



1.4 



[3 ~ 1.15. Therefore, for small oscillations {[3 < 1.15) Theorem 
vides the best range for blow up, correspondingly, for large oscillations 



{(3 > 1.15) Theorem 1.5 gives a better range. The curve "theoretical 
scattering" is the same as in Figure 9.2 and is given by Thm. 1.1, see 



(9.9) and values pj. The numerical threshold (dotted curve) is given in 
Table 9.8, values p~ and p'j^ . All p values in this graph are normalized 
Note that for small oscillations the numerical threshold 



by Vl + (3^. 

dotted curve coincides with the "theoretical scattering" curve. 
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